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A; .  *  stretching  stiffnesses 

*  J 

a,b  *  plate  dimensions  in  x  and  y  directions 

B- •  *  bending-stretching  coupling  stiffnesses 

*  J 

Crs  =  coefficients  defined  in  Eqs.  (2.2.4),  (3.3.4),  and  (4.3.4) 

Dy  =  bending  stiffnesses 

dx  =  3(  )/ax 

E  -  Young's  modulus  of  isotropic  ordinary  material 

Ec,Et  =  compressive  and  tensile  Young’s  modul i  (isotropic  bimodulus) 

k  k 

Eu,E22  ~  Young's  moduli  in  directions  x  and  y  (orthotropic  bimodulus) 

G13.G23  =  longitudinal-thickness  and  transverse-thickness  shear  moduli 

Gzc,G2£  -  transversely  isotropic,  bimodulus-material  shear  moduli 

' — b.. - —V.'  ,rV  total  thickness  of  plate 

i  V*'"* 

1^,..  rotatory  inertia  coefficient  per  unit  mid-plane  area 

2 

-  *4^X5  *  shear-correction  coefficients 

_ L  .  .  *  linear  differential  operators  defined  in  Eqs.  (2.1.10), 

--  (3.1.10),  and  (4.1.11) 

M^.Nj  *  stress  couples  and  inplane  stress  resultants 

m|,nJ  *  thermally  induced  stress  couples  and  inplane  stress 

,  1  resultants 

P  *  normal  inertia  coefficient  per  unit  mid-plane  area 

Qx.Qy  a  thickness-shear  stress  resultants 
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mid-plane  area 

■  temperature 

•  temperature  coefficients  defined  in  Eqs.  (4.4.1)  and  (4.4.2) 
«  time 
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■  displacements  in  x,y,z  directions 

■  mid-plane  displacements  in  x,y,z  directions 

•  fiber  volume  fraction 

■  matrix  volume  fraction 

■  bending-slope  coefficients  (amplitudes  of  ♦  ,♦  ) 

x  y 
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thickness  directions 

■  Vh-  Vh 

■  neutral -surface  positions  associated  with  e  »0  and  c  «0 

x  y 

■  coefficient  of  thermal  expansion 

•  ir/a,  ir/b 

•  strain  component  at  arbitrary  location  and  at  mid-plane 
«  Poisson's  ratio  of  Isotropic  material 

■  fiber  and  matrix  Poisson's  ratios 

«  major  (longitudinal-transverse)  and  transverse-thickness 
Poisson's  ratios 

•  stress  components 

■  slope  functions  in  x  and  y  directions 
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CHAPTER  I 


INTRODUCTION 

The  rate  of  progress  of  technological  innovation  is  dependent 
on  the  development  of  new  and  better  materials.  The  new  and  rapidly 
developing  composites  made  a  significant  impact  on  the  engineering  field 
and  are  responsible  for  the  tremendous  progress  that  has  been  achieved 
recently  in  the  structural  and  aerospace  industries. 

Composites  are  materials  made  up  of  more  than  one  constituent 
material.  According  to  this  literal  definition,  almost  all  materials 
used  in  civil  and  mechanical  engineering  are  composites.  Wood  consists 
of  lignin  and  cellulose  fibers  and  is  clearly  a  "natural"  composite,  but 
so  too  are  cast  iron,  steel  and  other  metallic  alloys,  brick,  natural 
stone,  and  of  course  reinforced  concrete.  The  fact  that  none  of  these 
materials  are  perfectly  isotropic  leads  us  to  closer  definition  of  a 
composite  as  understood  today,  especially  in  the  advanced  technology 
industries  such  as  aerospace  and  automotive. 

Composites  are  generally  laminates  in  which  a  matrix  material 
Is  reinforced  in  either  one  or  more  planes  with  filaments,  fibers  or 
fibrous  material,  giving  the  composite  enhanced  mechanical  properties 
over  those  of  either  the  matrix  or  the  reinforcement  when  used  alone. 

The  matrix  can  consist  of  metal,  ceramics,  glass,  concrete,  gypsum,  or 
resins,  and  the  reinforcement  can  be  metal  rods  or  filaments,  whiskers 


1 


of  silicon  carbide  or  nitride,  sapphire,  carbon  fiber,  boron  fiber,  and 
various  types  of  glass,  asbestos,  and  cellulose  fibers. 

Glass-reinforced  plastics  are  being  used  extensively  and 
successfully  in  the  manufacture  of  storage  vessels.  They  are  also  be¬ 
ginning  to  be  used  structurally  in  buildings,  and  a  good  deal  of  thought 
is  going  into  the  design  of  G.R.P.  bridges.  The  aerospace  industry  con¬ 
tinues  to  lead  in  the  use  of  composites  for  very  high-performance  appli¬ 
cations  with  products  ranging  from  rocket  casings  and  major  portions  of 
fuselage,  wing,  and  empennage  assemblies  to  compressor  blades  and  heli¬ 
copter  rotor  blades.  ' 

The  principal  reasons  for  using  composites  in  place  of  conven¬ 
tional  materials  are: 

(1)  Composites  are  anisotropic;  so  in  order  to  get  the  greatest 
economy  of  material,  either  for  cost  or  weight  saving,  the  reinforcing 
fibers  can  be  oriented  in  the  plane  where  they  will  be  most  effective. 

(2)  Composites,  unlike  metals,  can  often  be  molded  with  a  varying 
thickness  at  no  extra  cost.  This  gives  an  additional  freedom  to  econo¬ 
mize  the  material. 

(3)  Composites  have  improved  strength  and  stiffness,  especially 
when  compared  with  other  materials  on  a  unit  weight  basis.  For  example, 
composites  can  be  made  that  have  the  same  strength  and  stiffness  as 
high-strength  steel,  yet  are  70  percent  lighter!  Other  advanced  com¬ 
posites  are  as  much  as  three  times  as  strong  as  aluminum,  the  common 
aircraft  structural  material,  yet  weigh  only  60  percent  as  much! 

(4)  Composite  materials  can  be  tailored  to  efficiently  meet  design 
requirements  of  strength,  stiffness,  fatigue,  thermal  conductivity. 
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corrosion  resistance,  and  other  parameters  all  in  various  directions. 

The  advent  of  advanced  fiber-reinforced  composites  has  been 
called  the  biggest  technical  revolution  since  the  jet  engine.  This 
claim  is  very  striking  because  the  tremendous  impact  of  the  jet  engine 
on  military  aircraft  performance  is  readily  apparent.  The  impact  on 
commercial  aviation  is  even  more  striking  because  the  airlines  switched 
from  propeller-driven  planes  to  all -jet  fleets  within  the  span  of  just 
a  few  years. 

Currently,  almost  every  aerospace  company  ;  developing  pro¬ 
ducts  made  with  fiber- reinforced  composite  materials.  After  passing 
through  the  different  stages  of  usage,  people  are  now  dreaming  for  the 
final  stage  of  an  all-composite  high-performance  airplane. 

As  the  applications  of  fiber-reinforced  composites  in  structures 
become  more  widespread,  the  prediction  of  behavior  of  plates  constructed 
of  such  materials  become  increasingly  important.  One  of  the  character¬ 
istics  of  certain  composite  materials,  known  as  bimodulus  materials,  is 
that  they  exhibit  quite  different  elastic  properties  when  loaded  along 
the  fiber  direction  in  tension  as  opposed  to  compression  [1-4] (see  Fig.  1.1). 

These  materials  are  listed  in  Table  1.1.  In  the  literature, 
this  class  of  materials  has  variously  been  called  bilinear,  bimodulus, 
different-modulus,  and  multi -modulus.  Here  the  term  bimodulus  is  be¬ 
lieved  to  be  most  descriptive  of  a  material  having  different  linear 
stress-strain  relations  in  compression  than  in  tension. 

The  first  multi-dimensional  model  was  proposed  by  Ambartsumyan 
[  5]  for  Isotropic  material,  such  as  a  composite  material  with  spherical 
particles.  It  was  later  extended  to  the  orthotropic  case  [63- 
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The  second  and  third  models  are  the  restricted-compl iance 
model  due  to  Isabekyan  and  Khachatrayan  [7]  and  the  first-invariant 
model  of  Shapiro  [8].  A  fourth  model  is  the  weighted -compliance  theory 
originated  by  Jones  [9]. 

The  fifth  model  is  the  fiber-governed  bimodulus  symmetric 
compliance  model  originated  by  Bert  [10]. 


Fig.  1.1.  Bimodulus  idealization- 


A  plate  subjected  to  a  loading  which  produces  plate  bending  or 
vibration  obviously  experiences  both  tension  and  compression;  therefore, 
a  more  accurate  analysis  should  take  this  into  consideration. 
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Table  l.l  Some  Bimodulus  Materials 


Reinforcement 

Geometry 

Ref. 

Tensile  Young's  Modulus 
Divided  by 

Compressive  Younq's  Modulus 

ATJ-S  graphite 

Granular 

5 

1.2 

ZTA  graphite 

Granular 

5 

0.8 

Glass-epoxy 

Fibrous 

5 

1.25 

Boron-epoxy 

Fibrous 

5 

0.8 

Graphite-epoxy 

Fibrous 

5 

1.4 

Carbon-carbon 

Fibrous 

5 

2.0  to  5.0 

Kevlar-rubber 

Fibrous 

6 

0.77  (transverse)  * 
to  305  (longitudinal) 

Polyester-rubber 

Fibrous 

• 

0.75  (transverse)  * 
to  16.7  (longitudinal) 

Based  on  experimental  results  reported  by  Patel  et  al .  [2]. 


The  existing  literature  available  in  English  on  bending  of  bi¬ 
modulus  plates  is  quite  sparse  and,  with  only  a  few  exceptions,  is 
limited  to  bimodulus  isotropic  material  [8,11-14].  Shapiro  [8]  considered 
the  very  simple  problem  of  a  circular  plate  subjected  to  a  pure  radial 
bending  moment  at  its  edge,  but  he  used  Love's  stress-function  formula¬ 
tion  rather  than  plate,  theory.  Kamiya  [11]  treated  large  deflections 
(geometric  nonlinearity)  of  uniformly  loaded,  clamped-edge  circular 
plates,  using  an  iterative  finite-difference  technique.  In  [12],  Kamiya 
applied  the  energy  method  to  large  deflections  of  simply  supported  rec¬ 
tangular  plates  subjected  to  sinusoidally  distributed  loading.  In  [13], 
Kamiya  Included  the  effect  of  thickness  shear  deformation,  but  only  for 
the  simple  one-dimensional  case  of  cylindrical  bending.  The  only  analysis 
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applicable  to  anisotropic  bimodulus  material  is  the  work  of  Jones  and 
Morgan  [14],  who  treated  cylindrical  bending  of  a  thin,  cross-ply  lami¬ 
nate. 

it 

In  the  realm  of  plates  laminated  of  ordinary  anisotropic  mat¬ 
erials,  the  theory  due  to  Reissner  and  Stavsky  [15],  is  generally  recog¬ 
nized  as  the  classical,  linear  (small-deflection)  thin-plate  theory. 
Although  there  have  been  numerous  approximate  solutions  of  this  theory, 
only  a  relatively  few  closed-form  solutions  have  appeared.  Notable  among 
these  are  the  works  of  Whitney  [16]  and  Whitney  and  leissa  [17]  for  both 
antisymmetric  cross-ply  and  antisymmetric  angle-ply  rectangular  plates 
with  certain  (different)  kinds  of  simply  supported  edges.  For  an  in- 
infinitely  long  strip  of  finite  width,  Padovan  [18]  presented  a  solution 
for  the  case  of  an  arbitrary  laminate. 

Kamiya  [19,20]  considered  problems  of  thermal  stresses  in 
a  bimodulus  thin  plate.  An  annulus  with  axi symmetric  steady  temperature 
distribution  was  analyzed  numerically.  Ambartsumyan  [21,22]  presented 
a  general  theory  of  strains  and  stresses  for  bimodulus  materials  loca¬ 
ted  In  a  temperature  field.  Other  literature  available  on  thermal  bend¬ 
ing  [23-27]  deals  with  plates  of  ordinary  materials. 

Vibration  of  plates  has  been  treated  by  several  authors  [28- 
32],  but  the  problem  of  bimodulus  plate  vibration  has  not  been  attempted 
previously. 

Apparently,  the  present  work  is  first  to  consider  anisotropic, 

*  Throughout  this  report,  the  term  ordinary  will  be  used  to  distinguish 
materials  that  do  not  exhibit  bimodulus  action,  i.e.,  materials  in  which 
the  tensile  and  compressive  stiffnesses  coincide. 
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bimodulus,  thick  plates  finite  in  both  directions  in  closed  form  except 
a  few  exceptions  [33.343- 

The  problems  of  static  bending,  free  vibration,  thermal  bending 
and  non-linear  large  deflection  (of  thick  bimodulus  composite  rectangular 
plates)  have  been  analyzed  and  are  presented  separately  in  Chapters  II, 
III,  IV,  and  V.  Numerical  computations  have  been  carried  out  and  were 
compared,  and  good  agreement  was  obtained  with  existing  solutions  of 
special  cases  existing  in  the  literature. 
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CHAPTER  II 


STATIC  BENDING  OF  THICK  RECTANGULAR  PLATES 

Consider  the  case  of  ordinary  (not  bimodulus)  material.  A 
single-layer  plate  constructed  of  an  ordinary  material  that  is  macro- 
scoplcally  homogeneous  is  obviously  symmetric  about  the  midplane  of  the 
plate,  and  thus  there  is  no  coupling  between  bending  and  stretching 
during  small -deflection  bending.  Likewise,  a  plate  consisting  of  multi¬ 
ple  layers  of  ordinary  materials  of  various  thicknesses  arranged  sym¬ 
metrically  about  the  midplane  has  no  bending-stretching  coupling  at 
small  deflections.  However,  in  the  case  of  a  general  laminate,  i.e., 
one  not  symmetric  about  the  midplane,  bending-stretching  coupling  is 
induced. 

Now,  consider  the  case  of  a  single  layer  of  bimodulus  mat¬ 
erial.  The  different  properties  in  tension  and  compression  cause  a 
shift  In  the  neutral  surface  away  from  the  geometric  midplane,  and  sym¬ 
metry  about  the  midplane  no  longer  holds.  The  results  of  this  is  that 
a  single-layer  bimodulus-material  plate  exhibits  bending-stretching 
coupling  of  the  orthotropic  type  I.e.,  analogous  to  a  two-layer  cross- 
ply  plate  (one  layer  at  0°  and  the  other  at  90°)  of  ordinary  orthotropic 
material.  (See  Figs.  2.1,  2.2) 

Using  Bert's  fiber-governed  symmetric -matrix  macroscopic 
material  model  [10], it  can  be  assumed  that  there  are  two  symmetric  plane- 
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(c)  “wo-layer  cross-ply  witn  90“/0"  (a)  Two-layer  c-oss-ply  with  0°/90° 

stacking  sequence  stacking  sequence 

*1g.  2.1.  Laminate  configurations  ^or  rectangular  elates. 
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Jx  bottom 


(a)  Single  layer  with  1cngit~ain.il  fibers. 


Jy  bottom 


(b)  Single  layer  with  transverse  fibers. 


{ d)  Two-layer  cross-ply  laminate  with  0°  ?0°  stcokiv;  jeoience 


Pis.  l.  l.  Stress  uistr-outi  -  -s  l‘i  vdwlus  plaes. 
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stress  reduced  stiffness  matrices:  one,  when  the  fibers  are  in  tension 
along  their  length,  and  another,  when  they  are  in  compression  in  the 
same  direction.  Invoking  the  Voigt  hypothesis  in  the  fiber  direction, 
for  which  it  Is  well-established,  it  is  assumed  that  the  fiber-direction 
normal  strains  in  the  fibers  and  in  the  matrix  are  identical.  Then  the 
criterion  for  changing  from  tension  to  compression  can  be  taken  to  be 
the  fiber-direction  normal  strain  in  each  layer.  This  is  a  much  more 
convenient  criterion  to  apply  than  is  the  fiber-direction  normal-stress 
criterion. 


2.1  Governing  Equations 

Consider  a  plate  of  thickness  h  composed  of  an  even  number  of 
indentical  orthotropic  layers  bonded  together,  arranged  alternately  at 
angles  0°  and  90°.  The  origin  of  a  Cartesian  coordinate  system  is 
located  within  the  central  plane  (x-y)  with  the  z-axis  being  normal  to 
this  plane  (see  Fig.  2.3). 


Fig.  2.3.  Cartesian  coordinates  for  rectangular  plate. 


The  stress  resultants  and  stress  couples,  each  per  unit 
length,  are  defined  in  the  usual  way  as 


h/2 

(Nx*  V  Nxy’  V  V  *  l.h/2  (ox*  V  T*y’  Txz’  V)dz 


(M  ,  M  ,  M  ) 
'  x  y  xy' 


h/2 

<V  V  Txv)z  dz 

-h/2  x  y  xy 


(2.1.1) 

(2.1.2) 


The  theory  developed  by  Yang,  Norris,  and  Stavsky  [35]  is 
based  on  the  following  assumed  displacement  field, 

u  =  u°(x,y)  +  z<i>x(x,y) 

v  a  v°(x,y)  +  z*y(x,y)  (2.1.3) 

w  =  w(x,y) 


where  u,  v,  and  w  are  the  displacement  components  in  the  x,  y,  and  z 


and 


'  N  ' 

■ 

X 

N., 

y 

N 

J  xy 

►  = 

M 

X 

M, 

y 

M 

l  xyJ 

L 

directions,  respectively,  and  i|>  and  i|>  are  called  the  slope  functions. 

x  y 

The  constitutive  equations  for  an  unsymmetric  cross-ply  lam¬ 
inate  can  be  written  as  follows: 

An  A 1 2  0  Bn  B 1 2  0 

A12  A22  0  B12  B22  0 

0  0  A  66  0  0  B  66 

Bn  B 12  0  D11  D 1 2  0 

B 12  B22  0  D 1 2  D22  0 

0  0  B  66  0  0  D  66 


’y 


# 


x,x 


* 


y.y 


Uj  ^  lii 

\y,x  vx,y 


1. 


K4A44 


*5A55 


w,  ,  +  1 

< 

y  y  1 

/ 

l  w’x  +  *x  J 

(2.1.4) 


(2.1.5) 
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Here,  differentiation  is  denoted  by  a  comma,  and  the  exten- 
sional,  flexural -extensional  coupling,  and  flexural  or  twisting  stiff¬ 
nesses  for  a  laminate  of  an  arbitrary  number  of  layers  are  defined  by 

h/2 

“  \  ,z,z2)dz  ,  i  ,j*l  ,2,6  (2.1.6) 

-h/2 


In  addition  to  performing  the  integrations  in  a  piecewise 
manner  from  layer  to  layer,  one  also  has  to  take  into  account  the  pos¬ 
sibility  of  different  properties  (tension  or  compression)  within  a  layer. 
This  is  worked  out  in  detail  for  a  two-layer  cross-ply  laminate  in 
Appendix  A. 

2  2 

The  quantities  K4  and  K5  are  shear  correction  coefficients 
which  may  be  calculated  by  various  static  and  dynamic  methods  [36]. 

Taking  into  account  the  shear  deformation,  one  can  write  the 
equations  of  equilibrium  (neglecting  the  body  forces)  as  follows: 


N  +  N  *  0 

x,x  xy,y 

N  +  N  *  0 
xy.x  y,y 

Q  +  Q  *  -  q 
x,x  yy,y  y 

M  +  M  -0*0 

x,x  xy.y  vx 


M  +  M  -  Q  *  0 
xy.x  y,y  yy 


(2.1.7) 


Here  Nx>x  «  3Nx/ax,  N  *  aN^/ay,  etc. 

Substituting  equations  (2.1.4)  and  (2.1.5)  Into  equations 
(2.1.7),  we  obtain  the  equations  of  equilibrium  in  terms  of  the 
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generalized  displacements. 

A 1  lU  *xx  +  A66U*yy  +  (Ai2  +  Aee)v?xy  +  *u*XtXX  +  B66*x>yy 


+  (Bj2  +  B66)^y|Xy  *  ° 


(A12  +  A6e)u?xy  +  a66V?xx  +  A22v°yy  +  (B12  +  See)'!' 


x,xy 


+  Bs6^y,xx  +  B«*y,*y  =  0 


K5A55(*x>x  +  w,xx)  +  K24A44(^y>y  +  w,yy)  =  -  q 

(Bl2  +  B66)u?xy  +  B66v?xx  +  B22v°yy  +  (012  +  D66)^x>xy 


(2.1.8) 


+  D6«*y,xx  +  D22#yfyy  '  K4CA44Uy  +  W,y)]  =  0 


Bnu,xx  ♦  B66u?yy  ♦  (B12  ♦  B66)v?xy  ♦  Dn*x>xx  ♦  066*x>yy 
+  (D12  +  ■  K5CA55(lilx  +  w,x)3  =  0 


Or,  in  operator  form 

l-hu-J 


f  0 
u 

0 

< 

O 

0 

►  =  - 

q  [ 

% 

0 

hip 

k  x; 

°  J 

M=l,2,3,4,5 


(2.1.9) 


where  1s  a  metric  linear  differential  operator  matrix  with 
the  following  components: 
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Ill  =  Andx2  +  A66dy2 
L12  =  (A12  +  Agg )  dxdy 
L13  =  0 

Li4  =  C(Bi2  +  B6g)/h]  dxdy 
L15  s  (Bn/h)  d2  +  (B66/h)  dy2 
L22  =  Aggd^  +  A22dy2 
L23  =  0 

L24  ^  (Bgg/h)  dx2+  (B22/h)  dy2  (2.1.10) 

L25  =  Lm 

L33  =  -  K5A55dx2  -  K4A44dy2 
L34  2  -  Kj(A44/h)  dy 
L3S  5  "  Ks(A55/h)  dx 

L44  =  (Dgg/h2)  dx2+  (D22/h2)  dy2-  K4A44/h2 

L45  =  C(B12  +  Dgg)/h2]  dxdy 

L55  =  (Du/h2)  dx2+  (Dgg/h2)dy2  -  K^A55/h2 

2.2.  Application  to  Rectangular  Plate  Hinged  on  all  Edges 

The  boundary  conditions  for  a  rectangular  plate  simply  supported 
on  all  edges  can  be  specified  as  follows: 

Along  the  edges  at  x  *  0  and  x  =  a, 

W  !  ty  ’  Mx  =  0 

v°  -  Nx  *  0 
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Along  the  edges  at  y  *  0  and  y  3  b. 


(2.2.1) 


w  =  i>  3  M  *  0 
x  y 

u°  =  Ny  =  0 


Consider  the  loading  to  be  sinusoidally  distributed  in  both 
the  x  and  y  directions: 


Here 


q  *  q  sin  ax  sin  sy 
a  =  */a  ,  e  i  n/b 


Furthermore,  a  and  b  are  the  plate  dimensions  along  x  and  y  axes, 
respectively. 

The  governing  equations  (2.1.9)  and  boundary  conditions  (2.2.1) 
are  exactly  satisfied  in  closed  form  by  the  following  set  of  functions: 

u°  3  U  cos  ax  sin  ey 

v°  *  V  sin  ax  cos  ey 

w  *  W  sin  ax  sin  ey  (2.2.2) 

h^y  =  Y  sin  ax  cos  ey 

h^x  3  X  cos  ax  sin  ey 

Substituting  solutions  (2.2.2)  into  the  governing  equations 
(2.1.9),  one  obtains 


£Ck.^ 


u  ' 

r  0  ' 

V 

0 

w 

^  3  4 

V 

Y 

0 

X 

j 

1 

.  o. 

(2.2.3) 


k,«.3l  ,2, 3, 4, 5 
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where  [C  ]  is  a  5x5  symmetric  matrix  containing  the  following  elements: 

Ci i  =  -  Ana2  -  A66S2 
C i 2  =  -  (A12  +  A6g)a8 
Cl  3  =  0 

C i 4  =  •  [ ( B 1 2  +  B66)/^3a8 
C 1 5  i  -  (Bj i/h)a2  -  (B66/h)S2 
C22  =  "  Ag6a2  -  A22i32 

C23  =  0 

C24  s  -  (B66/h)a2  -  (B22/h)e2  ^.2A] 

C2S  5  C14 

C33  =  -  (K^Assa2  +  K^A4462) 

C34  =  ~  K4(A44/h)3 

C35  =  “  ^(Asg/hja 

C44  i  -  (D66/h2)a2  -  (D22/h2)32  -  K4(A44)/h2 

C4s  =  ~  C ( Dj 2  +  B5g)/hi]a8 

C55  =  -  (Du/h2)a2  -  (D66/h2)62  -  KsAss/h2 


2.3  The  Positions  of  the  Fiber-Direction  Neutral  Surfaces 
From  the  kinematics  of  the  deformation 


U,  +  Z  <|/ 

*x  x,x 


(2.3.1) 


V ,  +  Z  4* 

y  y.y 
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Thus,  the  neutral -surface  positions,  for  the  longitudinal  (x) 


and  transverse  (y)  directions,  respectively,  are 


hU/X 


-v° 

-  *  •  hv/Y 
y,y 


(2.3.2) 


So,  in  computations  the  values  for  znx  and  z  are  first 
assumed  to  obtain  the  displacements.  Actual  displacements  can  then  be 
obtained  by  an  iterative  procedure  with  the  help  of  the  equations 
(2.3.2). 


2.4  Numerical  Results 

As  the  first  example,  we  take  the  case  of  a  homogeneous 
(single-layer)  plate  of  transversely  isotropic  bimodulus  material.  The 
plane  of  isotropy  is  assumed  to  coincide  with  the  midplane  of  the  plate, 
and  the  inplane  Poisson's  ratio  is  assumed  to  be  zero.  Then  the  closed- 
form  solution  reduces  to  the  simplified  form  [33].  Numerical  results 
are  presented  in  Tables  2.1  and  2.2. 
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Table  2.1. 

Comparison  of  Neutral -Surface 
Transversely  Isotropic  Square 

Locations  for 
Plate 

Neutral  - 

Surface  Location  Zf 

Et/Ec’Gzt/G 

zc 

W0-1 

0.3 

0.5 

Exact  Closed-Form  Solution: 

0.5 

-  0.08578 

-  0.08578 

-  0.08578 

1.0 

0 

0 

0 

2.0 

+  0.08578 

+  0.08578 

+  0.08578 

Simplified  Approximate  Solution  [33] 

0.5 

-  0.08579 

-  0.08579 

-  0.08579 

1.0 

0 

0 

0 

2.0 

+  0.08579 

+  0.08579 

+  0.08579 

Mixed  Finite 

-Element  Solution  [33]: 

0.5 

-  0.08578 

-  0.08578 

-  0.08578 

1.0 

0 

0 

0 

2.0 

+  0.08578 

+  0.08578 

+  0.08578 

Table  2.2.  Comparison  of  Maximum  Deflections  for  Transversely 
Isotropic  Square  Plate  (b/h»0.l,  K- -b/G) 


Dimensionless  Deflection  WE.h‘/q  b** 

_  _  G  _  0 

G  /E  -0.1  0.3  0.5 

v  L  «.  L  »l  C 


Exact  Closed-Form  Solution: 


0.5 

0.05348 

0.04774 

0.04660 

1.0 

0.03688 

0.03283 

0.03201 

2.0 

0.02674 

0.02387 

0.02330 

Simplified  Approximate 

Solution  [33 

0.5 

0.05004 

0 . 046o0 

0.04591 

1.0 

0.03445 

0.03202 

0.03153 

2.0 

0.02530 

0.02342 

0.02296 

Mixed  F 1 n 1 

ite-Element  Solution  [33]- 

0.5 

0.05329 

0.04743 

0.04626 

1.0 

0.03675 

0.03261 

0.03178 

2.0 

0.02664 

0.02371 

0.02313 

It  is  noted  that  the  neutral-surface  location  Is  Independent  of 
GJc  and  G,t.  The  agreement  among  the  results  obtained  by  all  three  sol¬ 
utions  Is  quite  good. 

According  to  classical  thln-plate  theory  for  a  rectangular 
Isotropic  plate  with  simply  supported  edges,  the  dimensionless  maximum 
deflection  Is  given  by 


The  values  are  computed  for  three  aspect  ratios  using  the  above 
formula  and  are  compared  with  the  present  results  In  Table  2.3  below. 


20 


j 

j 


Table  2.3.  Dimensionless  Deflections  for  Rectangular  Isotropic 
Plate  as  Determined  by  the  Thin-Plate  Theory  and  the 
Present  Work 


WEi^/q^4 

Aspect  h/b»0. 1  h/b-0.01  h/b=0.001 

Ratio  - 


a/b 

Thick* 

Thin* 

Thick* 

Thin* 

Thick* 

Thin* 

0.5 

0.07439 

0.07392 

0.07401 

0.07392 

0.07392 

0.07392 

1.0 

0.02908 

0.02887 

0.02887 

0.02888 

0.02887 

0.02887 

2.0 

0.00478 

0.00462 

0.00462 

0.00462 

0.00462 

0.00462 

*  "Thick"  denotes  the  present  thick-plate  theory  and  "thin"  denotes 
classical  thin-plate  theory. 


As  examples  of  some  actual  bimodulus  materials,  aramid-cord/ 
rubber  and  polyester-cord/rubber  are  selected.  The  material  properties 
used  are  listed  in  Table  2.4.  The  data  are  based  on  test  results  of 
Patel  et  al. [2], using  the  data-reduction  procedure  of  [10],  except  for 
the  thickness  shear  moduli,  which  were  estimated. 
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Table  2.4.  Elastic  Properties  for  Two  Tire-Cord/Rubber, 
Unidirectional,  Bimodulus  Composite  Materials"1 


Property  and  Units 

Aramid- 

k-1 

■Rubber 

k-2 

Polyester-Rubber 

k-1  k-2 

Longitudinal  Young's  modulus,  GPa 

3.58 

0.0120 

0.617 

0.0369 

Transverse  Young's  modulus,  GPa 

0.00909 

0.0120 

0.00800 

0.0106 

Major  Poisson’s  ratio,  dimensionless4’ 

0.416 

0.205 

0.475 

0.185 

Longitudinal -transverse  shear  modulus,  GPa° 

0.00370 

0.00370 

0.00262 

0.00267 

Transverse-thickness  shear  modulus,  GPa 

0.00290 

0.00499 

0.00233 

0.00475 

a  Fiber-direction  tension  is  denoted  by  k«l,  and  fiber-direction  compression 
by  k-2. 

b  It  is  assumed  that  the  minor  Poisson's  ratio  is  given  by  the  reciprocal 
relation. 

j  It  Is  assumed  that  the  longitudinal-thickness  shear  modulus  is  equal  to 
this  one. 


Numerical  results  for  single-layer  rectangular  plates  with  the 
fibers  oriented  parallel  to  the  x  axis  are  given  in  Table  2.5,  while  those 
for  cross-ply  plates  are  listed  In  Table  2.6. 

The  exact  closed-form  solution  developed  here  for  thick,  rec¬ 
tangular  plates  of  single-layer  and  cross-ply  laminates  of  orthotropic  bi¬ 
modulus  material  has  been  shown  to  agree  well  with  an  existing  approximate 
solution  for  Isotropic  bimodulus  plates  and  with  a  mixed  finite-element 
solution. 


To  show  the  general  trend,  plots  have  been  presented  in  Fig.  2.4. 
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Table  2.5.  Neutral -Surface  Positions  and  Dimensionless  Deflections 
for  Rectangular  Plates  of  Single-Layer  0°  Aramid-Rubber 
and  Polyester-Rubber,  as  Determined  by  Two  Different 
Methods  (Thickness/Width,  h/b*0.1) 


Aspect 


WE"Ch3/q  b4 
_ 0 


Ratio 

C.F? 

F.E? 

C.F? 

F.E* 

C.F? 

F.E? 

Aramid-Rubber: 

0.5 

0.4453 

0.4454 

-  0.3304 

-  0.3007 

0.002544 

0.002750 

0.6 

0.4452 

0.4452 

-  0.2941 

-  0.2734 

0.004560 

0.004827 

0.7 

0.4447 

0.4447 

-  0.2564 

-  0.2419 

0.007393 

0.007712 

0.8 

0.4440 

0.4440 

-  0.2220 

-  0.2117 

0.01105 

0.01140 

0.9 

0.4431 

0.4431 

-  0.1923 

-  0.1846 

0.01545 

0.01582 

1.0 

0.4420 

0.4420 

-  0.1671 

-  0.1614 

0.02046 

0.02083 

1.2 

0.4394 

0.4394 

-  0.1285 

-  0.1250 

0.03160 

0.03193 

1.4 

0.4363 

0.4363 

-  0.1015 

-  0.09919 

0.04313 

0.04335 

1.6 

0.4328 

0.4329 

-  0.08228 

-  0.08070 

0.05406 

0.05416 

1.8 

0.4292 

0.4294 

-  0.06838 

-  0.06724 

0.06390 

0.06388 

2.0 

0.4253 

0.4254 

-  0.05813 

-  0.05727 

0.07250 

0.07236 

Polyester-Rubber: 

0.5 

0.3044 

0.3045 

-  0.1597 

-  0.1234 

0.001529 

0.001971 

0.6 

0.3044 

0.3045 

-  0.1538 

-  0.1245 

0.002652 

0.003265 

0.7 

0.3042 

0.3044 

-  0.1426 

-  0.1198 

0.004283 

0.005075 

0.8 

0.3039 

0.3041 

-  0.1299 

-  0.1124 

0.006517 

0.007487 

0.9 

0.3035 

0.3037 

-  0.1174 

-  0.1041 

0.009421 

0.01055 

1.0 

0.3029 

0.3031 

-  0.1061 

-  0.09586 

0.01303 

0.01430 

1.2 

0.3015 

0.3018 

-  0.08728 

-  0.08111 

0.02223 

0.02367 

1.4 

0.2999 

0.3001 

-  0.07329 

-  0.06941 

0.03348 

0.03492 

1.6 

0.2979 

0.2982 

-  0.06296 

-  0.06042 

0.04574 

0.04703 

1.8 

0.2957 

0.2960 

-  0.05528 

-  0.05356 

0.05793 

0.05897 

2.0 

0.2934 

0.2936 

-  0.04959 

-  0.04828 

0.06925 

0.07003 

*C.F.  denotes  closed-form  solution;  F.E.  signifies  finite-element  solution  [33]. 
(For  in-plane  displacements,  see  Appendix  D) 
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Table  2.6.  Neutral -Surface  Positions  and  Dimensionless  Deflections 
for  Rectangular  Plates  of  Two-Layer  Cross-Ply  Aramid- 
Rubber  and  Polyester-Rubber,  as  Determined  by  Two  Different 
Methods  (Thickness/Width,  h/b»0.1) 


Aspect 

Z 

X 

Zy 

WE2; 

'h3/qob4 

Ratio 

C.F? 

F.E? 

C?F? 

F.E? 

C.F? 

F.E? 

Aramid-Rubber: 

0.5 

0.4433 

0.4431 

-  0.06343 

- 

0.06223 

0.002472 

0.002576 

0.6 

0.4427 

0.4426 

-  0.05478 

- 

0.05443 

0.004388 

0.004518 

0.7 

0.4418 

0.4418 

-  0.04794 

- 

0.04778 

0.007072 

0.007220 

0.8 

0.4407 

0.4407 

-  0.04247 

- 

0.04237 

0.01054 

0.01070 

0.9 

0.4396 

0.4396 

-  0.03803 

- 

0.03795 

0.01475 

0.01490 

1.0 

0.4384 

0.4384 

-  0.03437 

- 

0.03430 

0.01957 

0.01972 

1.2 

0.4356 

0.4356 

-  0.02883 

- 

0.02860 

0.03043 

0.03054 

1.4 

0.4326 

0.4325 

-  0.02470 

- 

0.02477 

0.04185 

0.04190 

1.6 

0.4292 

0.4292 

-  0.02160 

- 

0.02165 

0.05282 

0.05280 

1.8 

0.4257 

0.4256 

-  0,01922 

- 

0.01923 

0.06277 

0.06264 

2.0 

0.4219 

0.4219 

-  0.01735 

- 

0.01734 

0.07151 

0.07137 

Polyester-Rubber : 

0.5 

0.3650 

0.3652 

-  0.1285 

- 

0.1256 

0.002539 

0.002732 

0.6 

0.3644 

0.3646 

-  0.1178 

- 

0.1164 

0.004527 

0.004772 

0.7 

0.3638 

0.3639 

-  0.1097 

- 

0.1089 

0.007288 

0.007575 

0.8 

0.3631 

0.3631 

-  0.1036 

- 

0.1031 

0.01078 

0.01109 

0.9 

0.3622 

0.3622 

-  0.09886 

- 

0.09859 

0.01487 

0.01519 

1.0 

0.3613 

0.3613 

-  0.09526 

- 

0.09502 

0.01933 

0.01966 

1.2 

0.3593 

0.3593 

-  0.09001 

- 

0.09000 

0.02846 

0.02879 

1.4 

0.3571 

0.3570 

-  0.08660 

- 

0.08660 

0.03674 

0.03707 

1.6 

0.3546 

0.3545 

-  0.08430 

- 

0.08430 

0.04356 

0.04389 

1.8 

0.3519 

0.3518 

-  0.08267 

- 

0.08267 

0.04890 

0.04925 

2.0 

0.3491 

0.3490 

-  0.08150 

. 

0.08150 

0.05301 

0.05337 

*C.F.  denotes  closed-form  solution;  F.E.  signifies  finite-element  solution  [33]. 
(See  Appendix  D  for  in-plane  displacements.) 


CHAPTER  III 

VIBRATION  OF  THICK  CROSS-PLY  LAMINATED 
BIMODULUS  RECTANGULAR  PLATES 

Thin-plate  theory  does  not  take  into  account  either  the  effect 
of  transverse  shear  deformation  or  rotatory  inertia,  and  hence  it  becomes 
inaccurate  for  thicker  plates.  Mindlin  [37]  considered  both  of  these 
effects  for  homogeneous  isotropic  plates,  by  assuming  that  the  displace¬ 
ment  variation  across  the  thickness  is  linear  for  u  and  v  and  constant 
for  w.  He  also  had  to  assign  a  value  to  the  shearing  rigidity  factor  on 
suitable  physical  considerations.  His  solution  does  not  satisfy  the 
governing  elasticity  equations  exactly,  but  does  permit  the  satisfac¬ 
tion  of  a  set  of  three  boundary  conditions  on  each  edge.  Mindlin, 
Schacknow,  and  Deresiewicz  [38]  applied  this  method  to  the  vibrations 
of  thick  rectangular  plates  with  two  opposite  sides  simply  supported  and 
the  other  two  edges  with  various  conditions.  The  present  work,  to  the 
author's  knowledge.  Is  the  first  to  consider  the  bi modulus  property  In 
the  thick  cross-ply  rectangular  plates. 

3.1  Governing  Equations 

Consider  a  plate  of  thickness  h  composed  of  an  even  number  of 
identical  orthotropic  layers  bonded  together,  arranged  alternately  at 
angles  0°  and  90°  (see  Fig.  2.2). 
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The  stress  and  moment  resultants,  each  per  unit  length,  are 
given  in  the  usual  manner  as 


r  h/2 

,N  ,Q  ,Q  )  * 

(a„ia„, dz 

(3.1.1) 

xy  x  y 

-h/2 

x  y  xy  xz  yz 

rh/2 

(M  ,M  ,M  )  ■ 

'  x  y  xy 

-h/2 

{Vay*Txy)z  dz 

(3.1.2) 

The  displacement  components,  u,  v,  and  w  in  the  x,  y,  and  z 
directions  respectively  can  be  expressed  in  terms  of  mid-plane  displace¬ 
ments  u°,  v°,  w°  and  slope  functions  and  as: 

x  y 

u  =  u°(x,y,t)  +  z<»x(x,y,t) 

v  *  v°(x,y,t)  +  Zif/y(x,y,t)  (3.1.3) 

w  *  w(x,y,t) 


where  t  is  time. 

Constitutive  equations  for  an  unsymmetric  cross-ply  laminate, 
as  has  already  been  mentioned  in  Chapter  II,  are: 


N. 


N. 


xy 


xy 


An 

Au 

0 

Bn 

Bl2 

0 

A1 2 

a22 

0 

B 1 2 

b22 

0 

0 

0 

A66 

0 

0 

Bee 

Bn 

B 1 2 

0 

Du 

Dl2 

0 

Bl2 

b22 

0 

°12 

D22 

0 

0 

0 

®66 

0 

0 

Dee 

0 

v,y 
v?x  + 


x,x 


yy.y 

\b  +  iii 

y.x  x,y ) 


and 


K4A44 


K5A5  5 


w,y  +  *y 

w’x  +  *x 


(3.1.4) 


(3.1.5) 
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Differentiation  here  Is  denoted  by  a  comma,  i.e.,  (  )»x  = 

3(  )/ax,  and  the  extensional,  flexural-extensional  coupling,  and  flex¬ 

ural  stiffnesses  for  the  laminate  are  defined  by 

h/2 

(V8irV  *  \.w  (3.1.6) 

1  ,j*l ,2,6 

As  usual,  K*  and  k|  are  shear-correction  coefficients. 

Taking  Into  account  the  shear  deformation  and  the  rotatory 
Inertia,  the  equations  of  motion  can  be  written  as  follows: 


Nx,x  +  Nxy,y  3  P  u,tt  +  R  ^x.tt 
Nxy,x  +  Ny,y  *  P  v,tt  +  R  *y,tt 
Qx.x  *  «y,y  ’  P  “•« 

Mx.x  *  Mxy,y  '  1«  '  Ru°tt  +  1  ♦x.tt 


*  N„  „  -  Q„  ■  «  *?,.  +  '  * 


xy.x  y.y  y 


tt  "'y.tt 


(3.1.7) 


Here  P,  R,  and  I  are  the  normal,  in-plane,  and  rotatory 

Inertia  coefficients  per  unit  mid-plane  area  and  are  defined  by 

h/2 

(P,R,I)  ■  p(l,z,z2)  dz  (3.1.8) 

'-h/2 


where  p  is  the  material  density. 

Substituting  equations  (3.1.4)  and  (3.1.5)  into  equations 
(3.1.7),  we  obtain  the  equations  of  motion.  In  operator  form. 
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f“0' 

0  j 

< 

0 

0  ! 

w 

s  ,  , 

0  V 

h*y 

0 

>x( 

0 

where  [l^J  Is  a  symmetric  linear  differential  operator  matrix  with  the 
following  elements: 

La  =  Andx2  +  A^2-  Pdt2 

L12  -  (A12  *  6  )  dxdy 

La  =  0 

LlH  -  t(®12  +  Bgg ) /h]  dxdy 

L1S  s  (Bn/h)dx2  ♦  (B66/h)dy2  -  (R/h)  d‘ 

L22  =  A^d  ;  *  A22dy2  -  Pd  2 

L23  =  o 

L24  =  (B66/h)dx2+  (B22/h)dy2  -  (R/h)  d‘  (3.1.10) 

L25  5  Llu 

L33  =  -  K25A55dx2  -  K2A44dy2  +  Pd2 
L34  -  ■  Ki»(A44/h)  dy 

L35  5  "  ^(Ajj/h)  dx 

L44  s  (0s6/h2)dx2+  (D22/h2)dy2-  K4A44/h2  -  (I/h2)dt* 

L„5  s  C(0l2  ♦  066)/h2]  dxdy 

L55  5  [Dn/h2)dx2+  (06b/h2)dy2-  K2AS5/h2  -  (I/h‘)dti 
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3.2  Application  to  Plate  Simply  Supported  on  all  Edges 
The  boundary  conditions  are: 


Along  the  edges  at  x  3  0  and  x  3  a, 

w  3  *y  =  Mx  3  0 

v°  3  N  3  0 
x 

Along  the  edges  at  y  3  0  and  y  3  b, 

w  =  *x  =  My=  0 

u°  3  N  =0 

y 


(3.2.1) 


3.3  Closed-Form  Solution 

The  governing  equations  (3.1.9)  and  the  boundary  conditions 
(3.2.1)  are  exactly  satisfied  in  closed  form  by  the  following  set  of 
functions: 

o  .  iwt 

u  3  U  cos  ax  sin  Sy  e 

v°  =  V  sin  ax  cos  sy  eiwt 

w  =  W  sin  ax  sin  sy  e1wt  (3.3.1) 

hij)y  3  Y  sin  ax  cos  Sy  e1wt 

hii»x  3  X  cos  ax  sin  sy  e1wt 

Here,  <*>  is  the  natural  frequency  associated  with  the  mode  having  axial 
and  transverse  wave  numbers  m  and  n,  and 

a  =  mtr/a  ,  S  =  mr/b  (3.3.2) 

where  a  and  b  are  plate  dimensions  in  the  x  and  y  directions,  respectively. 

Substituting  solutions  (3.3.1)  into  the  governing  equations 
(3.1.9),  we  obtain  the  following: 
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The  frequency  w  can  be  determined  by  setting  |C^J  *  0. 

3.4  Neutral -Surface  Locations 

This  is  the  same  as  in  the  preceding  chapter:  znx  =  -  hU/X; 
zny  *  "  ''terat''ve  procedure  is  used  to  obtain  the  final  dis¬ 

placements. 


3.5  Numerical  Results 

Since  there  is  no  previous  analysis  for  vibration  of  bimodulus 
plates,  the  present  one  could  be  compared  only  with  rectangular  plates 
laminated  of  ordinary  materials. 

Comparisons  with  Jones, and  Fortier  and  Rossettos  are  presented 
in  Tables  3.1  and  3.2  below.  It  can  be  seen  that  the  agreement  is  good. 


Table  3.1.  Comparison  of  Fundamental  Natural  Frequencies 
of  Rectangular  Antisymmetric  Cross-Ply  Plates 
at  Different  Plate  Aspect  Ratios 

(Eu/E22=40  ,  Gl2/E:2=0.5  ,  vu«0.25  ,  b/h=10) 


Aspect 
Ratio 
a/b  • 

ob:  r 

IT-  / 

Jones  [31] 

Present 

0.5 

2.050 

1.934 

1.0 

0.825 

0.794 

1.5 

0.650 

0.612 

2.0 

0.580 

0.565 

2.5 

0.560 

0.548 

3.0 

0.550 

0.541 
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Table  3.2.  Comparison  of  Fundamental  Frequency  (m=n"l) 
of  a  Square  Cross-Ply  Plate 

(E  \  i / E_. ^ =40  ,  G|  ./E?:=l ,  G23/E„-0.5.  ^*i^=0.25) 


a/h 

ua  </o/  E ,  2  h  ^ 

Fortier  and  Rossettos  [28] 

Present 

10 

10.80 

10.80 

50 

11.65 

11.65 

Typical  results  for  (bimodulus)  aramid-rubber  and  polyester- 
rubber  are  tabulated  in  the  following  tables  (Tables  3.3  and  3.4).  See 
Table  2.4  for  the  elastic  properties  and  Appendix  B  for  the  densities. 

Computations,  based  on  the  closed-form  solution,  have  been 
carried  out  for  thick,  rectangular  plates  of  cross-ply  laminates  and 
compared  with  existing  works.  Close  agreement  was  reached. 

Abrupt  changes  in  the  values  of  Zx  and  are  noticeable  which 
may  be  due  to  the  bimodulus  effect  in  combination  with  the  eigenvalue 
nature  of  this  problem  (see  Fig.  3.1). 


Table  3.3.  Values  of  Dimensionless  Fundamental  Frequencies  and 

Neutral -Surface  Locations  for  Single-Layer  Orthotropic 
Rectangular  Plates  Having  b/h=10  and  Different  Aspect 
Ratios 


a/b 

z  z 

x  y 

u>b^/p/D22/ A2 

Aramid-Rubber: 

0.4 

0.0100  0.0137 

14.2370 

0.6 

-  0.0204  0.0296 

8.9952 

0.8 

-  0.0346  0.0515 

6.3712 

1.0 

0.0818  -  0.3619 

4.1732 

1.2 

-  0.0724  -  0.2206 

3.8780 

1.4 

0.3145  0.4442 

2.9399 

1.6 

0.0058  -  0.5673 

2.4122 

1.8 

0.0179  -  0.3207  x 

10"4  2.1262 

2.0 

-  0.0257  -  0.0243 

1.8935 

Polyester-Rubber: 

0.4 

0.0679  -  0.2661 

12.8460 

0.6 

0.7121  x  10'4  0.9147  x 

10'5  7.4975 

0.8 

0.1314  0.0952 

4.8319 

1.0 

0.0471  -  0.0197 

3.3832 

1.2 

-  0.0613  0.0182 

2.5949 

1.4 

0.0103  0.0034 

2.1992 

1.6 

0.0903  -  0.1354 

3.2867 

1.8 

-  0.0178  0.0474 

1.9318 

2.0 

-  0.0178  0.0474 

1.4641 
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Table  3.4.  Values  of  Dimensionless  Fundamental  Frequencies  and 
Neutral -Surface  Locations  for  Cross-Ply  Rectangular 
Plates  Having  b/h-10  and  Different  Aspect  Ratios 


a/b 

Zx 

z 

y 

uib‘-rip/{)  jj/ir- 

Aramld- 

■Rubber: 

0.4 

-  0.2196  x 

10'2 

-  0.0152 

14.1740 

0.6 

0.1929  x 

10‘2 

-  0.0363 

8.5910 

0.8 

-  0.0205 

-  0.0353 

6.0253 

1.0 

-  0.0196 

-  0.0383 

4.4521 

1.2 

0.0312 

-  0.0209 

2.0092 

1.4 

-  0.0109 

0.0159 

2.6748 

1.6 

0.0142 

-  0.7792  x 

10'4  2.3851 

1.8 

-  0.1409  x 

10'2 

0.1801  x 

10"2  0.8723 

2.0 

-  0.0241 

0.3046  x 

10‘3  1.7419 

Polyester-Rubber: 

0.4 

0.0742 

0.4955  x 

10*2  7.7991 

0.6 

0.0283 

0.1847  x 

10'2  4.3010 

0.8 

-  0.0276 

0.0219 

2.7333 

1.0 

-  0.1028 

0.0449 

1.9735 

1.2 

-  0.1989 

0.0667 

1 . 6037 

1.4 

-  0.8590  x 

10’3 

0.1042 

1.2888 

1.6 

-  0.2123  x 

10'4 

0.3120 

1.133 

1.8 

-  0.3218  x 

10'5 

0.7538 

1.002 

2.0 

1.2780 

0.8525 

0.9613 
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Fig.  3.1.  Variation  of  fundamental  vibration  frequency  with  aspect 
ratio  for  two-layer  cross-ply  rectangular  plates. 


CHAPTER  IV 

THERMAL  BENDING  OF  CROSS-PLY  THICK 
BIMODULUS  RECTANGULAR  PLATES 

Based  upon  the  mathematical  theory  of  elasticity  of  bimodulus 
materials  and  upon  the  Neumann  hypothesis,  Ambartsumyan  [22]  dealt  with 
the  development  of  the  theory  of  thermoelasticity  for  elastic  bimodulus 
material.  Kamiya  [20]  developed  the  fundamental  equations  for  axisym- 
metric  plane  stress  problems  for  a  bimodulus  thin  plate.  Das  and  Rath  [26] 
presented  an  analysis  of  thermal  bending  for  a  moderately  thick  rectangular 
plate  subjected  to  a  temperature  distribution  which  is  antisymmetric  about 
the  middle  plane  of  the  plate,  but  is  arbitrary  along  the  direction  perpen¬ 
dicular  to  simply  supported  edges  and  constant  along  the  other  perpen¬ 
dicular  direction.  Bapu  Rao  [39]  treated  thermal  bending  of  thick  iso¬ 
tropic  rectangular  plates  taking  into  consideration  the  shear  deformation 
capability. 

The  present  analysis  deals  with  the  thermal  bending  of  unsym- 
metrically  cross-plied  bimodulus  rectangular  plates  simply  supported  on 
all  edges. 


4.1  Governing  Equations 

The  plate  Is  subjected  to  a  sinusoidal  temperature  distribution, 
T,  Thus,  temperature  terms  appear  in  the  constitutive  equations.  The 


37 


thermoelastic  constitutive  relations  for  the  material  are  written  as 
follows: 


r 

°x 

"Qn 

Ql2 

0 

(s) 

f 

«l  -  axT 

1  0y 

►  = 

Ql2 

Q22 

0 

€2  "* 

l  °xy. 

0 

0 

Q66_ 

v 

(4.1.1) 


The  stiffness  matrix  [Q]  takes  different  values  in  tension  and  compres¬ 
sion  depending  upon  the  sign  of  the  fiber-direction  strain. 


< 

\ 


ijc 

if  ei 

<  0 

ijt 

if  £i 

>  0 

(4.1.2) 


The  coefficients  of  thermal  expansion  <34  and  a2  in  the  x  and  y 
directions,  respectively,  also  depend  upon  the  sign  of  the  fiber-direction 
strain.  Also, 

f  air  lf  £i  <  0 

0  1  (4.1.3) 

1  ajt  if  ei  >  0 


The  laminate  constitutive  relations  can  be  represented  by: 


(  T 

N  +  N 

X  X 

Au  A12  0 

Bn/h 

B12/h 

0 

f  0  ^ 

u’x 

V"y 

^12  ^22  0 

Bu/h 

B22/h 

0 

V? 

y 

0  0  a66 

0 

0 

B6$/h 

V?  +u? 
x  ’y 

T 

>  ■ 

> 

(Mx  +  Mx)/h 

Bu/h  B12/h  0 

Dn/h2  Du/h2 

0 

hip 

*x,x 

Bu/h  B22/h  0 

012/h2  D22/h2 

0 

hip 

y*y 

.V”  . 

0  0  B66/h 

0 

0 

W2 

-J 

v  +  h*x  y 

(4.1.4) 
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where  the  thermally  induced  inplane  forces  are 


f 

i 


Qll“l  +  Ql2a2 

Ql2al  +  Q22a2 


T  dz 


(4.1.5) 


(4.1.6) 


(4.1.7) 


As  usual  the  stretching,  stretching-bending,  and  bending  stiff 
nesses  for  the  laminate  are  defined  as 


h/2 


•VW  ’  f  <V('-2-22) 

-h/2 


(4.1.8) 


i  ,j=l ,2,6 


Including  shear  deformation,  the  equations  of  equilibrium  are: 


1  +  N 

=  0 

X.x  xy,y 

1  +  N 

*  0 

xy,x  y,y 

1  +0  » 

■  0 

X.x  7  ,y 

1  +  M 

-  Q 

-  0 

x.x  xy,y 

yx 

1  +  M 

-  Q 

»  0 

xy.x  y.y 

y 

(4.1.9) 


Substituting  equations  (4.1.4)  and  (4.1.5)  Into  equations 
(4.1.9),  we  get: 
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r  t 

Nx,x 

0 

V 

Ny.y 

j  0 

►  =  4 

o  f 

MI,x 

h*v 

v  *  J 

mt 

(4.1.10) 


k,i*l ,2, 3, 4, 5 


where  [l^]  is  a  symmetrix  linear  differential  operator  matrix  with  the 
following  elements: 


Lll  5  AHdx2  +  As6dy 
Li2  -  (Ai2  +  A66)  dxdy 


Li  3  =  0 

Li4  5  (Bi2  +  B66)  dxdy 
Li 5  5  Bndx  +  Begd^ 

L22  =  A66dx  +  A22d^ 


i-23  -  0 

L24  3  B^d2  +  B^d2  (4.1.11) 

1-25  =  Ll4 

L33  2  -  K^gd2  -  K^d2 
2 

L34  3  “  Ki,A44d^ 

L35  =  -  KsAssdx 

L44  =  066dx2+  ^22dy2“  K4A44 

L45  3  (0l2  +  ^66)  dxdy 

L55  s  Dndx2+  D66dy2-  KsA55 
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4.2  Application  to  a  Simply  Supported  Plate 
The  boundary  conditions,  as  usual,  are: 

Along  the  edges  at  x  ■  0  and  x  =  a, 

w  “  V  Mx  +  =  0 

v°  *  Nx  +  Nx  =  0 

Along  the  edges  at  y  *  0  and  y  *  b,  (4.2.1) 

w  =  iii  3  M  +M^  =  0 

x  y  y 

u°  “  Ny  +  Ny  =  0 

4.3  Closed-Form  Solution 

The  governing  equations  and  the  boundary  conditions  can  be 
satisfied  exactly  by  the  following  for  T  sinusoidally  distributed  along  x  and  y: 

u°  *  U  cos  ax  sin  By 

v°  3  V  sin  ax  cos  ey 

w  3  W  sin  ax  sin  ey  (4.3.1) 

h<j)y  3  Y  sin  ax  cos  ey 

hij>x  3  X  cos  ax  con  ey 

where 

a  3  mir/a  ,  e  3  nir/b  (4.3.2) 

Here  m  and  n  are  Integers,  and  a  and  b  are  plate  dimensions  in  the  x  and 
y  directions. 

Inserting  solutions  (4.3.1)  into  the  governing  equations 
(4.1.10),  we  get  the  following: 
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'ktJ 


u  ' 

f  r 
"I.X 

V 

Ny,y 

<  w 

>  *  < 

0  r 

y 

«;.x 

,x  , 

.  My.y , 

(4.3.3) 


where 

Cn  - 
C 1 2  = 
cu  = 
Cu  = 

C15  2 

C2;  = 

C2  3  = 
C24 
C;  5 

C33 
C34 

C35 


k,tal ,2, 3, 4, 5 

[Cj^]  Is  a  symmetric  matrix  containing  the  following  elements: 

-  An**2-  A666- 
“  (Aj2  +  A66)a6 
0 

“  (Bj2  +  B66 )at3 

-  Bna2  -  B6662 

*  Agga4,  -  A226‘ 

0 

-  Bgga‘  -  B226-  (4.3.4) 

Cu 

-  kjAjja-  -  K^Ai^S^ 

-  K4A448 


KsAssa 


C44  s  *  Ogga5  -  02282  -  K4A44 
C45  5  -  (012  +■  0gg)a8 

C$S  5  "  Ollai  "  °66fl2  “  Ks^ss 
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4.4  Mean  Temperature  and  Temperature  Gradient  Sinusoidally 
Distributed  over  a  Rectangular  Region 


where 


T(x,y,z)  ■  TQ(x,y)  *  (z/h)  T^x.y) 


TQ(x,y)  «  Tq  sin  (nwx/a)  sin  (mry/b) 
T^x.y)  *  sin  (mux/a)  sin  (n*y/b) 


(4.4.1) 


(4.4.2) 


For  Case  I, 


z  >  0  and  z  <  0  with  z  governing 
nx  ny  nx 


layer  1  (0°)  and  zny  layer  2  (90°)  (see  Appendix  C  for  the  remaining  cases). 


T  ny 

<;  ■  [  (q 


1122  a122  +  Q1222  a222)  TdZ 


U 

(Q1112  a112  +  Q1212  a212}  TdZ 


+  J  (Q1121  a121  +  Q1 221  a221)  TdZ 
0 

h/2 

+  ,  (Q1111  “ill  +  Q1211  a211)  TdZ 


(Q1122  a122  +  Q1222  °222)  *  8122 
(qill2  a112  +  qi212  a212*  “  ell2 
(qll21  a121  +  ql221  a221 ^  *  8121 
(Qllll  “ill  +  Q12U  “211*  *  8111 


(4.4.3) 


Then, 


NI  *  S122  To(zny  *  h/2>  *  s112  To<°  -  V  +  sl21  Vz„x  '  °> 

*  S1 1 1  To(h'2  '  z„x>  *  s122<V2h)(zJy  -  "2/4> 

*  8„2(T,/2h)(0  -  zzy)  *  ai21(Tl/2h)(z=x  -  0) 

*  8,„<T,/2h)(h=/4  -  zzx) 

Nx  ’  *s122  *  B1 1 1  > <T0h/2)  *  (a121  -  8,,,)  T0znx  +  (S,22  -  8n2) 

ToZny  +  ^8111  ‘  3i22^Tlh/8)  +  ^ S1 21  '  *1 11  ^TlZnx/2h* 

+  U122  -  3112)(T12^y/2h)  (4.4.4) 

Similarly, 

Ny  =  (t3222  +  B211)(T0h/2)  +  (6221  -  B211)  TQZnx  +  (B222  -  B212)  TQZny 

+  (6211  -  B222)(T1h/8)  +  (g221  -  B211)(TlZnx/2h)  +  -  6212) 

(Vny/2h)  {4-4-5) 

Now, 

h/2  0  z  h/2 

T  i-  nx  r  ' 

M  «  81?7  Tzdz  +  6-1,0  Tzdz  +  3,o,  Tzdz  +  3,,,  Tzdz 

J.h/2  J2ny  H2  J0  121  j2nx  111 

*  (8111  ”  8122^oh2/8)  +  (8121  '8in^oznx/2^  +  ^ S1 22  '  8112^ozny/2^ 

+  (S122+ 8ul)(T1h2/24)+(8121  -  6ul)CT1zJx/3h)+(6122- 

(4.4.6) 
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And  similarly. 


My  *  (S211  '6222)fToh2/8)*(6221  '  S211 }  Vnx/2) +(8222  "  Vny/2) 

+  (s222  +  62n  )(r1h2/24)+(e221  -  e2n  )^ri2nx/3h)+  ( e222 ' 82  Wny^ 

(4.4.7) 

Using  equations  (4.4.4),  (4.4.5),  (4.4.6),  and  (4.4.7),  we  obtain  the 

following: 

■I.x  *  a{(8122  +  8lll^oh/2)  +  (8121  ‘  81 1  Voznx  +  (s122  '  6112) Vny 
+  (0m  “  e122)^lh/8)  +^121  ■  6lll)<flznx/2h) 

♦  (0122  -  0n2^flzny/2h)}  {4’4-8) 

Ny,y  *  6{(e222  +  a211^oh/2)  +  (e221  "  S211 }  Vnx  +  (e222  “  S212)  Vny 
+  (0211  -  0222^1h/8)  +  (6221  '  82lVVnx/2h) 

+  ^222  "  8212^1zny^2^  (4.4.9) 

Mx,x  *  a{(8lll  '  8122^oh2/8)  +  (8121  '  8lll^Vnx/2)  +  ^122  "  8112) 

(Vny/2^  +  (6i22  +  8lll^^lh2/24^  +  (e121  “  8111  ^^lznx/3h) 

+  (0122  •  8112)(Vny/3h)}  (4.4.10) 

and 
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S{(32n  -  e222)(Toh~/8)  +  (6221  -  B211^Vnx/2)  +  (S2 22  '  6212^Vny/2) 
+  (e222  +  B21 1 ) (^1h“/24)  +  (*221  ‘  8211^Vnx/3h) 

+  (g222  -  S212)(T1zny/3h)  (4.4.11) 


4.5  Neutral-Surface  Locations 
As  explained  already, 

znx  =  -  hU/X  ,  zny  =  -  hV/Y  (4.5.1) 

Numerically,  the  values  z  and  z  are  computed  as  follows: 

*  n  Jr 

The  values  of  z  and  z  are  assumed  in  the  beginning  to  get 
ma  ny 

displacements.  The  above  equations  are  used  to  obtain  the  new  values  of 

z„„  and  z„  and  fed  back  to  get  one  a  more  accurate  set.  This  procedure 
nx  ny 

is  repeated  until  the  actual  deflections  are  obtained. 


4.6  Numerical  Results 

Numerical  results  are  compared  with  Boley  and  Weiner's  work  [40] 
for  an  isotropic-ordinary-material,  single-layer,  thin  rectangular  plate. 
Close  agreement  has  been  obtained  (see  Table  4.1  below). 


Table  4.1.  Comparison  with  Boley  and  Weiner's  Work  [40]  for  an  Isotropic 
Single-Layer  Thin  Rectangular  Plate  at  Different  Aspect  Ratios 
(Ejj/E22  =  1.00  ,  v^2  =  V21  3  0.3  ,  b/h  3  10) 


Aspect 

Ratio, 

a/b 

Deflection,  W/h 

(m=n=l ) 

8oley  and  Weiner  [40] 

Present 

0.5 

0.5300 

0.5264 

1.0 

6.5858 

6.5789 

1.5 

6.3112 

6.3063 

2.0 

2.1104 

2.1058 
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Typical  numerical  results  for  (bimodulus)  aramld-rubber  and 
polyester-rubber  are  listed  In  the  following  tables  (see  Table  2.4  for 
the  properties). 

The  solution  for  an  isotropic, ordinary-material, single-layer, 
thin  rectangular  plate  has  been  specialized  from  the  present  analysis  and 
Is  compared  with  such  a  solution  available  In  the  literature.  Good  agree 
ment  was  obtained. 

Sudden  change  In  the  deflection  has  been  observed  in  the  case 
of  aramld-rubber.  To  show  the  general  trend,  graphs  have  been  plotted 
(see  Fig.  4.1). 
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Table  4.2.  Values  of  Wh/a*  ^b2  and  Neutral -Surface  Locations 
for  a  Single-Layer  Orthotropic  Rectangular  Plate 
[b/h“10  ,  aiV<MC*0.5  »  oi2^/a2C=l  .0  ,  aiVoi2^“Q.l  , 

t0/T 1*1 .0) 


a/b 

Zx 

zy 

Wh/a^  Tjb2 

0.5 

-  2.3007 

Aramid-Rubber: 

-  0.1666 

0.4533  x  10'1 
0.5547  x  10"1 

0.6 

-  1.8521 

-  0.2237 

0.7 

-  1.5573 

-  0.2745 

0.6579  x  10'1 

0.8 

-  1.3482 

-  0.3160 

0.7636  x  10'1 

0.9 

-  1.1192 

-  0.3473 

0.8722  x  10"1 

1 .0 

-  1.0708 

-  0.3690 

0.9839  x  10"1 

1.2 

-  0.8959 

-  0.3886 

0.1218 

1.4 

-  0.7766 

-  0.3861 

0.1467 

1.6 

-  0.6914 

-  0.3709 

0.1732 

1.8 

-  0.6286 

-  0.3498 

0.2011 

2.0 

-  0.5819 

-  0.3266 

0.2301 

0.5 

-  2.1099 

Polyester-Rubber: 

-  0.6511 

0.1401  x  10'1 

0.6 

-  1.6790 

-  0.8482 

0.1764  x  10'1 

0.7 

-  1.4029 

-  1.0080 

0.2156  x  10*1 

0.8 

-  1.2112 

-  1.1218 

0.2583  x  10'1 

0.9 

-  1.0710 

-  1.1903 

0.3046  x  10'1 

1.0 

-  0.9649 

-  1.2194 

0.3546  x  10_1 

1.2 

-  0.8172 

-  1.1933 

0.4672  x  10'1 

1.4 

-  0.7233 

-  1.1062 

0.5940  x  10'1 

1.6 

-  0.6630 

-  1.0015 

0.7318  x  10_1 

1.8 

-  0.6258 

-  0.9016 

0.8757  x  10'1 

2.0 

-  0.6060 

-  0.8158 

0.1020 

(See  Appendix  0  for  the  in-plane  displacements) 
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Table  4.3.  Values  of  Wh/a2t  T^2  and  Neutral -Surface  Locations 

for  a  Cross-Ply  Rectangular  Plate  (b/h*10  ,  aifya-^O.l  , 
aiVaiC*0.5  *  a2^/a2C“l*0  »  T  /f j “  1 . 0 ) 


a/b 

Zx 

.  z* 

Wh/a^  T^2 

Aramld-Rubber: 

0.5 

0.1276 

-  0.1510  x  103 

-  0.2620  x  10'1 

0.6 

0.1184 

-  0.1722  x  103 

-  0.3983  x  10’1 

0.7 

0.1165 

-  0.2347  x  103 

-  0.5973  x  10'1 

0.8 

0.1191 

-  0.5316  x  103 

-  0.8814  x  10'1 

0.9 

0.0578 

-  0.3450 

0.3660  x  10'2 

1.0 

0.0980 

-  0.3498 

0.4432  x  10'2 

1.2 

0.1304 

-  0.4327 

0.5848  x  10‘2 

1.4 

0.1399 

-  0.3375 

0.6998  x  10‘2 

1.6 

0.1402 

-  0.3339 

0.7852  x  10“2 

1.8 

0.1358 

-  0.3316 

0.8453  x  10’2 

2.0 

0.1286 

-  0.3302 

0.8862  x  10'2 

Polyester-Rubber: 

0.5 

0.7442 

-  0.8415 

-  0.5453  x  10'1 

0.6 

0.7433 

-  0.8470 

-  0.7863  x  10'1 

0.7 

0.7477 

-  0.8516 

-  0.9269  x  10'1 

0.8 

0.7584 

-  0.8439 

-  0.8915  x  10'1 

0.9 

0.7772 

-  0.8237 

-  0.6903  x  10*1 

1.0 

0.8075 

-  0.7973 

-  0.3825  x  10'1 

1.2 

0.9274 

-  0.7628 

0.2514  x  10'1 

1.4 

3.2603 

0.4369 

-  0.7672 

1.6 

3.3521 

0.4656 

-  0.8295 

1.8 

3.3178 

0.4636 

-  0.8473 

2.0 

3.2126 

0.4440 

-  0.8402 

(See  Appendix 

0  for  the 

In-plane  displacements) 
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(a)  Kev Tar-Rubber 


(b)  Polyester-Rubber 

Variation  of  dimensionless  deflection  with  respect  to 
aspect  ratio  for  two-layer  cross-ply  rectangular  plates 


CHAPTER  V 


LARGE  DEFLECTIONS  OF  BIMODULUS  CROSS-PLY 
THIN  RECTANGULAR  PLATE 

This  class  of  problem  has  been  attempted  In  different  ways  by 
various  people  [11,12,41-49].  Nonlinear  bending  of  ordinary  orthotropic, 
single-layer  thin  plates  subjected  to  uniform  loading  has  been  treated  by 
Niyogi  [41]  among  others. 

Perhaps  the  first  large  deflection  analysis  of  unsymmetrlcally 
laminated  plates  was  due  to  Plster  and  Dong  [42],  who  considered  isotropic 
rectangular  plates.  The  arbitrarily  laminated  fully  anisotropic  equations 
of  the  von  Karman  type  were  probably  first  presented  by  Whitney  and  Leissa 
[17],  who  did  not  solve  them.  Large-deflection  analyses  of  unsymmetrical , 
laminated  rectangular  plates  have  been  published  in  [43-49]. 

Large-deflection  analyses  of  plates  made  of  bimodulus  materials 
have  been  limited  to  isotropic  bimodulus  materials.  Kamlya  [11,12]  ana¬ 
lyzed  both  the  circular  and  the  rectangular  planforms. 

Due  to  the  complicated  algebra  Involved  In  the  nonlinear  be¬ 
havior  of  thick  plates,  the  work  in  this  chapter  has  been  reduced  to  thin 
plates.  An  approximate  solution  Is  obtained  by  using  the  Galerkln 
technique. 
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5.1  Basic  Equations 

Consider  a  thin  rectangular  plate  of  thickness  h  and  plate 
dimensions  a  and  b  (in  the  x  and  y  directions,  respectively)  subjected 
to  nonlinear  bending. 

In  view  of  Kirchhoff's  thin-plate  hypothesis,  the  displacement 
components  u,  v,andw  in  the  x,  y,  and  z  directions  can  be  expressed  in 
terms  of  mid-plane  displacements  u°,  v°,  and  w°  as: 

u  *  u°(x,y)  -  z  w,x(x,y) 

v  *  v°(x,y)  -  z  w,y(x,y)  (5.1.1) 

w  =  w(x,y) 


where  the  comma  denotes  differentiation. 

The  laminate  constitutive  relations  can  be  written  as: 

( 


"y 

1  N 

xy 

Mx 

My 

."xy 

Au  A12  0 


B11  B12  0 


aJ2  a22  o  b12  b22 


0  a66  o  o  b 


66 


Bn  b12 


'll  u12 


B12  B22  ^  B12  B22  0 


0  0  b66  0  0  0 


66 


“•x  + 


v?,  +  W  ,2/2 


U°  + V?  +  w,  w, 
’y  ’x  ’x  y 


^  (5.1.2) 


w, 


XX 


■  *’yy 

-  2w, 


xy 


The  quantities  A^,  B^,  and  0^  are  the  stretching,  stretching- 
bending  coupling,  and  bending  stiffnesses  defined  by 
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(A« 


(5.1.3) 


h/2 

'1j,B1j,D1j*  *  |_h/2  ^*z'z2)  Qij  dz 

1  ,j*l ,2,6 

The  equilibrium  equations  are  (neglecting  body  forces  and  body  moments): 


N  +  N  *0 
x,x  xy,y 

N  +  N  ■  0 
xy,x  y,y 

Nw,  +  2N  w,  +  N  w,  +  M  +2M  +M  +  q  «  0 
x  ’xx  xy  ’xy  y  *yy  x,xx  xy.xy  y,yy  M 


(5.1.4) 


where  q  represents  the  normal  load. 

Substituting  equations  (5.1.2)  Into  equations  (5.1.4),  we  get 
the  following.  For  equilibrium  in  the  x  direction: 

All(«?xx  +  W-*W’yJ  +  A12(v°*v  +  W‘vW’*V>  '  BHW»yxx  * 


’X  XX' 


xy 


y  xy' 


xyy 


+  ^ °'yy  +  u,xy  +  w’xw*yy  +  "VV  *  2B**  w’xyy  ‘  0 


or 


Anu?xx  +  ft66u?yy  ♦  (Au  +  Agjv?^  +  w,x(Anw.xx  +  A66w,yy) 

♦  w,y(A12  +  A66)w,xy  -  Bnw,xxx  -  (B12  +  2B66)w.xyy  »  0  (5.1.5) 

For  equilibrium  In  the  y  direction: 

A“(u°^  *  v°xx  *  w>x"’xy  *■  "•y">xx)  '  2B«6"-xxy 

♦  *  ‘w^y,  *  >*•/%,> 

‘  8‘*"'«xy  '  8^"'yyy  *  0 


or 
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(5.1.6) 


A22V,yy  +  A66v ‘XX  +  (A12  +  A66)u?xy  +  W,y(A22W,yy  +  AggW,^ ) 

*  w,x(A12+A66)W,xy  -  e22w,yyy  -  (B12  +  2B66)w,_  =  0 
For  equilibrium  in  the  2  direction: 


xxy 


[An(u?x  +  JW.2)  +  A12 (v?y  t  I„,2)  t  8U(.  »,xx)  *  Bu(-  w. 


’yy^w,xx 


+  2CA66(u?y  +  v?x  +  w,xw,y)  +  Bgg  ( -  2w ,¥U)>, 


xy  ’xy 


*  4  «•«)  *  |  %2>  -  B.2%x  -  S22“vy>-y, 


+  Bu<U’xxx  *  "-x"’xxx  *  »4>  +  6I2(,exxy  *  w.yw.xxj  *  w.^1 


"  °11W,XXXX  “  °12W> 


xxyy 


+  2[B66(u’xyy  +  V’xxy  +  W’xW’xyy  +  w’xy  +  w’y*Sxy  +  w-xx«>y) 

"  2D66w,xxyrI  +  B12(u°xyy  +  w’xw»xyy  +  w.Jy^ 

+  B22(v,yyy  +  w,yw,yyy  +  w,yy}  '  °i2w>xxyy  -  °22%yyy  +  q  =  0 


or 


BllU*XXX  +  (B12  +  2B66)(u?xyy  +  V?xxy)  +  B22v? 


yyy 


+  W»XCBUW*XXX  +  fB12  +  2B66)wiXyy]  +  W,v[B22  W,vvv  +  (B12  +  2B66)w,vvu] 


yyy 


xxy- 


+  2%x(B66  -  B12)w,  +  2*»*xyw»xV(B12  +  B66)  -  w,  B22w, 


yy  xy  *xy 

■  Dliw’xxxx  -  2(°12  +  2D66)w*xxyy  '  °22W*yyyy 

+  {<,  +  W‘Xx[An(uV?W‘i)  +Ai2(v?y  4w,2y)] 

lw  2, 

2  WV 


yy  “  yy 


•y  T  2  w’y 

+  w,yy[Ai2(u?x  +  1  w.2)  +  A22(v°y  +  I  w,J)] 

+  2w’XytA66(u»y  +  V«x  +  w,xw*y)  *  2B66w»vu^  *  0 


xy" 


(5.1.7) 
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5.2  Simply-Supported  Boundary  Conditions 

Along  the  edges  at  x  *  0,4 

w  ■  w.  *  M  *  0 
*yy  x 

v°  »  Nx  -  0 

Along  the  edges  at  y  ■  0,b  (5.2.1) 

w  '  "-XX  ■  My  ■  0 

u°  *  Ny  *  0 


5.3  Solution 

The  set  of  equations  (5.1 .5-5.1 .7)  are  coupled  and  nonlinear  in 
nature,  and  an  exact  solution  appears  to  be  extremely  difficult  to  obtain. 
Hence,  an  approximate  solution  will  be  obtained  here.  Let 
q  «  Q  sin  ax  sin  Sy 
w  »  W  sin  ax  sin  sy 

0  (5.3.1) 

u  ■  U  cos  ax  sin  sy 

v°  *  V  sin  ax  cos  sy 

By  the  Galerkin  method  (equations  (5.1 . 5) -( 5. 1 .7)  and  equations  (5.3.1)  are 
combined) : 

|  |  {(AnU»XX  +  A6*v?yy  +  (A12  +  A66)v?^  +  W,x(AUW,xx  ♦  A66w,yy) 

♦w,y(Ai2  ■•‘A66)>»,xy -Bnw,xxx  -  (B12  +  2866)w,  )>  cos  ax  sin  sy  dxdy-0 

or 
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a  b 

|  |  (’  An  p  u  cos2  aX  sin2  -  A66  p  U  cos2  ax  sin2  6y 
0  0 

-  V(A12  +  A66)  p  cos2  ax  sin2  6y  -  An  p-  W2  cos2  ax  sin3  sy  sin  ax 


or 


or 


or 


3 

-  A66  jgz  W2  cos2  ax  sin3  gy  sin  ax 

_3  3 

+  W2(A12  +  A66)  jgy  cos2  ax  cos2  8y  sin  ax  sin  ey  +  BnW  p  cos2  ax  sin2  ey 
+  W(B12  +  2B66)  cos2  ax  sin2  ey}  dxdy  ■  0 


-  An  Jr  u  x  “  a66  57  u  x  "  ^Al2  +  A&6*v  ab  T 

-  An  JT  w2*fp  -  a66  w2’f£  +  w2(Al2  +  Aee)  ib* 

+  BnWp-f  +  W(B12  +  2B66)  0 


U(+ 


Anir2b  A66ir2a 


-)  +  VO  (A12  +  A66) 


4a  4b 

Bnff3b  w3 

-  wt+  “35 —  +  (Bu  +  2B66) 

+  w2[+  |  AX1  Jy  +  |  a66  J  |  ( a i 2  +  a66)  j~]  *  o 


U{&  (Aub2  +  A66«2)>  +  V{x  (A12  +  A66)}+W{-  0  Bu-g  (Bl2  +  2B66) } 


+  W2  { |  ^  (Au  p  +  A66  -  2  (A12  +  A66))}  »  0 


b2 


(5.3.2) 
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Similarly, 


(A22v?yy  +  A66v°xx  +  (AU  ♦  A66)u°xy  +  *,y(H22*'yy  ♦  A66w,xx) 

♦  w,v(A12  +  A66N,  -B22w,  -(B12  +  2B66)w,  >sin  ax  cos  ey  dxdy-0 


a  u 

*  2  2 

{-  A22  p  V  sin2  ax  cos2  ey  -  A66  p  V  sin2  ax  cos2  ey 

0  0 

TT2  3 

■  (fti2  +  A66)  U  ^  sin2  ax  cos2  ey  -  A22W2  p  sin3  ax  cos2  ey  sin  ey 

IT3 

-  A66W2  pg-  sin3  ax  cos2  ey  sin  ey 

♦  (A12  +  A66)W2  sin  ax  cos2  ax  sin  ey  cos2  ey  +  B22  p  W  sin2  ax  cos2  ey 
+  (Bl2  +  2B66)  W  sin2  ax  cos2  ey}  *  0 


-*»&*¥ -  *„£v  u£$-AMu*££ff 


it3  8ab 


’  a66«2  ffr  +  ( A 1 2  +  a66)w2  |p  +  b22w  p  f- 


jt3  ab 


♦  (Bu  ♦  2B66)W  £  f  *  0 


v(+  A22  ^  +  A66  $)  +  U(Al2  +  A66)  £ 

-  W[B22  $+  (8U  ♦  2B66)  ji] 

+  W2[+  A22  fgl  +  A66  -  (A12  +  A66)  »  0 
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or 


u<V  (A12  +  A66)}  +  V{i-  (§a22  +  f  a66)} 
+  Ht-  B22  0  -  (B12  +  2B66)} 


+  W^9a  ^22  +  A&6  "  \  (Ai2  +  a6s)}  3 


(5.3.3) 


Also, 


a  b 

t 


\  \  {BllU*XXX  +  (B12  +  2B66)(u?  +  V°  )  +  B22v? 

0  0  ■ 


yyy 


+  w’x[Bnw’XXx  +  <Bi2  +  2B66)w,  ] 


xyy 


+  w.vCB22w.yyy  +  (B12  +  2B66)w,vvJ 


xxyJ 


+  2W*XX(B66  -  Bl2)w.yy  +  2W*XVW’XV  {B12  +  B66)  -  W>VVB22W, 


xy  ’xy 


yy  “  yy 


'  D11W,XXXX  ■  2(D12  +  2D6e)W»xxyy  -  022w’ 


yyyy 


+  q  +  W*XX[Au(U°x  +  l  w'^  +  A12(v?u  +  l  w *)] 


’X  2  ’x 

1  ..  2 ' 


y  2  y 
i  ..  2, 


+  w*yy  tAl2(«.x  +  2  w’x>  +  A22(v°y  +  2  W’y^ 

+  2w’xy  tA66(u»y  +  v?x  +  w,xw,y ^  "  2B66w>xy]>  sin  <*x  sin  sy  dxdy=0 


or 
a  b 


/»  3  3 

J  tjr  Bn  U  sin2  ax  sin2  8y  +  (B12  +  2B66)(^t  U  +  ^  V)  sin2  ax  sin*  By 
0  0 


+  B22  p-  V  sin2  ax  sin2  sy  -  Bnw2  pr  cos2  ax  sin  ax  sin3  sy 
-  (B12  +  2B66)  W2  cos2  ax  sin  ax  sin3  sy 
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"  B22W2  pr  sin3  <**  sin  By  cos2  sy 

„4 

-  (B12  +  2B6e)  W2  sin3  ax  sin  sy  cos2  By 

4 

+  2(B66  -  B12)  pp-W2sin3  ax  sin3  By 

4 

+  2W2{B12  +  B66)  sin  <*X  sin  sy  cos2  aX  cos2  By 

4  4 

-  B22  pr  W2  sin3  ax  sin3  ey  -  0UW  V  sin2  ax  sin2  sy 

4  4 

-  2(012  +  2D66)W  sin2  ax  sin2  sy  -  D22W  sin2  ax  sin2  By 
+  Q  sin2  ax  sin2  By  +  Au  73-  UW  sin3  ax  sin3  By 

a 

4  3 

-  An  2^tr  W3  sin2  ax  cos2  ax  sin4  By  +  A12WV  sin3  ax  sin3  By 

4  3 

-  A] 2  W3  sin4  ax  cos2  By  sin2  By  +  A12  UW  sin3  ax  sin3  By 

IT4  IT3 

■  a12  2d2b7  W3  cos2  aX  S3°2  aX  sin4  6y  +  Az2  b7  VW  si°3  aX  sin3  ^ 

4 

-  A22  W3  sin4  ax  cos2  By  sin2  By 

n3 

+  2A66  UW  sin  ax  cos2  ax  sin  sy  cos2  By 

IT3 

+  2A66  VW  sin  ax  cos2  ax  sin  By  cos2  By 

4 

+  2A66  W3  sin2  ax  cos2  ax  sin2  By  cos2  By 

4 

-  4B66  W2  ^7^7-  sin  ax  cos2  ax  sin  By  cos2  By>  *  0 
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nl  a  if  sb  /_  \/^3U  tt3Vv  t3  ab  „  u2  it4  8ab 

p-  tsHu  4  +  (B12  +  2B66)(-^-  +  -^-)  +  B22  V  —  -  BnWz  ^  jct 


HT  V  T  -  «„»  p  ^: 


o  u2  Z1*  8a b  / d  j.  id  \  7,14  i,7  8a b  n  tt4  8a b 

BuW  p  gp  (812  +  2B66)  pzW2^-  B22W-  p  ^p 


-  (b12  +  2B66)  pp  W2  gp  +  2W2(B66  -  B12)  pp 


tt4  16ab 


-  2W2(B12  +  B66)  pp  |p  -  B22  pp-  W2  ^  -  DnW  p  ^ 

-  2(D12  +  2D66)  W  ^  ^  -  D22W  j£^+Q^  +  An£uW^ 


n  n4  3ab  .  „  7t 3  16ab  „  it4  3ab  ,  „  it3  ,,,,  16at 

-  An  2p  W  -gp  +  A12WV  pg  -g p-  -  A12  2pp  W3  -p-  +  A12  pr  UW  -gp 


ir3  16ab 


-  fli2“32?rp‘^  +  *22  ET  ™  -  *22  55T“3^+  2*66  UW  jp-f 

*  2*66™  “3  3^§  -  4B«  “2  pps?'11 

6(6,,  4jT  +  (Bl2  +  2B66  ifo)  +  V((8,2  +  2B&6 )  ^  +  B2.  ^4-} 

*«-  »110-  2(0l2  +  2°66)  ^5  -  !>22  0) 

*  «2l-Bu  -  (B,2  *  2B66)  g  -  B22  ^  -  (SU  ♦  2B66) 

*  2(B66  -  S12)  ^  *  2(6,2  6  B65)  -  b22  ^  -  4B66  |^) 

*  w3<‘*“  TOSa*  '  *12  TS8a¥  "  *12  T28aE  ■  *22  mb1  *  2*66  rfab1 


it3  4ab 


.  I  h  1/  A  1  Sb^T  ,  m  1  6-TT  m  Qtt  1 

+  UW{An  -gp-  +  A12  -gg-  +  A66  g^} 

2.  tt  1 6  .  «  ir!6a  ,  w4,  .  „  ab  _  n 

+  VW{A12  -gp  +  A22  -gp-  +  2A66  9^>  +  %  X  '  0 


(5.3.4) 


The  above  three  equations  yield  a  cubic  equation  in  W  which  can 
be  solved  by  using  a  standard  subroutine  (such  as  ZRPOLY). 
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5.4  Neutral -Surface  Positions 


By  Kirchhoff's  hypothesis,  the  neutral-surface  locations  are 

defined  by  e  =  e  =0  which  implies 
J  x  y 


znx  3  uVw*xx  3  U/Wa 


Zny  =  Vw'yy  =  v/we 


(5.4.1) 


The  iterative  method  is  applied  as  in  the  other  problaits  to  get  the  actual 
deflection. 


5.5  Numerical  Results 

The  present  solution  is  compared  with  that  of  Kamiya  [11],  the 
only  solution  available  for  large  deflections  of  bimodulus  Isotropic 
rectangular  thin  plates, simply  supported  on  all  edges.  Good  agreement  is 
obtained  (see  Table  5.1). 


Table  5.1.  Comparison  of  Nondimensional  Deflection,  W/h,  with 
Kamiya's  Solution  [11].  (a/b=l,  vc*0.2,  q  aVE;^4 

=16.91 ) 


Et/Ec 

16.91 

WE22ChVqoa“ 

Z  *z  /h  Z  =z  /h 
x  nx  y  ny' 

Kamiya 

Present 

Present 

1.0 

0.4000 

0.4010 

-0.0795  -0.0795 

1.5 

0.3200 

0.3177 

-0.0398  -0.0100 

2.0 

0.2700 

0.2658 

0.2658  -0.0167 
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Table  5.2.  Dimensionless  Deflections  and  Neutral -Surface  Locations 
for  a  Single-Layer  Orthotropic  Rectangular  Plate  at 
(i)  qQ  -  qobVE22Cb4  3  1.0  ,  (ii)  qQ  -  qob4/E22ch4  =  200 


a/b 

2 

X 

(D  ... 

1 

X 

(ii) 

- 1 

(i) 

I - 

(ii) 

W/h 

(i) 

W/h 

(ii) 

Aramid- 

-Rubber: 

0.5 

0.0768 

0.0876 

-  0.0688 

-  0.0803 

0.6651 

X 

10"4 

0.0139 

0.6 

0.0780 

0.0828 

-  0.0636 

-  0.0724 

0.1384 

X 

10*3 

0.0278 

0.7 

0.0780 

0.0732 

-  0.0598 

-  0.0667 

0.2567 

X 

10'3 

0.0489 

0.8 

0.0795 

0.0611 

-  0.0564 

-  0.0644 

0.4377 

X 

10"3 

0.0774 

0.9 

0.0795 

0.0462 

-  0.0536 

-  0.0659 

0.7001 

X 

10~3 

0.1126 

1.0 

0.0805 

0.0290 

-  0.0512 

-  0.0699 

0.1063 

X 

10~2 

0.1530 

1.2 

0.0805 

-0.0093 

-  0.0477 

-  0.0885 

0.2178 

X 

10'2 

0.2435 

1.4 

0.0805 

-0.0498 

-  0.0464 

-  0.1146 

0.3945 

X 

10‘2 

0.3395 

1.6 

0.0805 

-0.0906 

-  0.0446 

-  0.1447 

0.6528 

X 

10’2 

0.4365 

1.8 

0.0789 

-0.1370 

-  0.0446 

-  0.1765 

0.1003 

X 

10’1 

0.5326 

2.0 

0.0773 

-0.1699 

-  0.0446 

-  0.2086 

0.1449 

X 

10'1 

0.6270 

Polyester-Rubber: 

0.5 

0.0706 

0.0565 

-  0.0515 

-  0.0441 

0.3090 

X 

10'3 

0.0572 

0.6 

0.0706 

0.0376 

-  0.0491 

-  0.0390 

0.6361 

X 

10'3 

0.1058 

0.7 

0.0706 

0.0126 

-  0.0471 

-  0.0390 

0.1166 

X 

10'2 

0.1674 

0.8 

0.0716 

0.0126 

-  0.0454 

-  0.0504 

0.1962 

X 

10'2 

0.2389 

0.9 

0.0716 

-0.0447 

-  0.0441 

-  0.0619 

0.3087 

X 

10'2 

0.3094 

1.0 

0.0716 

-0.0747 

-  0.0441 

-  0.0796 

0.4597 

X 

10‘2 

0.3832 

1.2 

0.0703 

-0.1334 

-  0.0427 

-  0.1231 

0.8974 

X 

10'2 

0.5305 

1.4 

0.0683 

-0.1899 

-  0.0427 

-  0.1712 

0.1529 

X 

10’1 

0.6743 

1.6 

0.0656 

-0.2433 

-  0.0437 

-  0.2202 

0.2344 

X 

10_1 

0.8130 

1.8 

0.0623 

-0.2929 

-  0.0472 

-  0.2681 

0.3298 

X 

10'1 

0.9453 

2.0 

0.0589 

-0.3382 

-  0.0503 

-  0.3137 

0.4337 

X 

10'1 

1.0700 

(See  Appendix  D  for  the  in-plane  displacements) 
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Table  5.3.  Dimensionless  Deflections  and  Neutral -Surface  Locations 
for  a  Cross-Ply  Bimodulus  Rectangular  Plate  at 
(1)5  =  qrtb4/E22ch1+  =  1.0  .  (11)  q  «  q  bVE22ch4  =  200 


a/b 

X 

.  (D  _ 

X 

(ii) 

(i) 

i.l) 

W/h 

_ _ (i) 

W/h 

(11) -  -  _ 

Aramid- 

■Rubber: 

0.5 

0.0934 

0.0875 

-  0.0280 

-  0.0316 

0.7149 

X 

10‘4 

0.0139 

0.6 

0.0934 

0.0827 

-  0.0247 

-  0.0329 

0.1480 

X 

10'3 

0.0278 

0.7 

0.0934 

0.0731 

-  0.0221 

-  0.0398 

0.2734 

X 

10-3 

0.0488 

0.8 

0.0934 

0.0611 

-  0.0202 

-  0.0544 

0.4645 

X 

io-3 

0.0774 

0.9 

0.0934 

0.0462 

-  0.0189 

-  0.0766 

0.7402 

X 

10'3 

0.1126 

1.0 

0.0934 

0.0292 

-  0.0189 

-  0.1033 

0.1121 

X 

IO*2 

0.1531 

1.2 

0.0934 

-0.0072 

-  0.0172 

-  0.1633 

0.2284 

X 

10"2 

0.2404 

1.4 

0.0921 

-0.0397 

-  0.0172 

-  0.2184 

0.4127 

X 

10‘2 

0.3187 

1.6 

0.0909 

-0.0625 

-  0.0172 

-  0.2580 

0.6800 

X 

10‘2 

0.3741 

1.8 

0.0895 

-0.0752 

-  0.0172 

-  0.2815 

0.1041 

X 

10*1 

0.4077 

2.0 

0.0874 

-0.0826 

-  0.0188 

-  0.2950 

0.1498 

X 

10'1 

0.4270 

Polyester-Rubber: 

0.5 

0.0895 

0.0634 

-  0.0815 

-  0.0959 

0.3519 

X 

10'3 

0.0616 

0.6 

0.0895 

0.0438 

-  0.0805 

-  0.1128 

0.7206 

X 

10‘3 

0.1124 

0.7 

0.0895 

0.0181 

-  0.0805 

-  0.1413 

0.1312 

X 

IO*2 

0.1759 

0.8 

0.0895 

-0.0104 

-  0.0805 

-  0.1789 

0.2188 

X 

10'2 

0.2469 

0.9 

0.0895 

-0.0396 

-  0.0805 

-  0.2216 

0.3403 

X 

io'2 

0.3307 

1.0 

0.0895 

-0.0678 

-  0.0805 

-  0.2673 

0.4999 

X 

10'2 

0.3931 

1.2 

0.0869 

-0.1437 

-  0.0832 

-  0.4042 

0.1527 

X 

IO'1 

0.5994 

1.4 

0.0869 

-0.1437 

-  0.0832 

-  0.4042 

0.1527 

X 

IO'1 

0.5994 

1.6 

0.0854 

-0.1554 

-  0.0856 

-  0.4356 

0.2217 

X 

IO'1 

0.6444 

1.8 

0.0840 

-0.1576 

-  0.0882 

-  0.4517 

0.2947 

X 

IO*1 

0.6669 

2.0 

0.0828 

-0.1542 

-  0.0909 

-  0.4605 

0.3659 

X 

IO'1 

0.6776 

(For  the  in-plane  displacements,  seee  Appendix  D) 


Comparison  of  the  present  theory  with  the  existing  solution  for 
bimodulus  isotropic  thin  plates  shows  close  agreement.  Typical  computa¬ 
tions  are  shown  in  the  above  two  tables,  and  graphs  are  presented  (Fig. 
5.1)  to  observe  the  general  trend. 
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0  TOO  200  300  400 
q  =  q0a4/E22V 

(a)  Isotropic  Et/Ec  =  2,  vc  *  0.2 


q  -  qoa4/E22ch4 

(b)  Aramid-Rubber 
(Cross-Ply) 


q  »  q0a4/E22Ch4 

(c)  Polyester-Rubber 
(Cross-Ply) 


Fig.  5.2.  Variation  of  dlmensionaless  deflection  with  dimensionless 
load  for  bimodulus  square  plate  (a/b*l).  (Dashed  line 
represents  linear  case.) 
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CHAPTER  VI 


CONCLUSIONS 

Good  agreement  exists  between  the  closed-form,  small-deflection 
solutions  presented  here  and  previous  approximate  analyses  carried  out 
by  several  authors  for  special  cases  (isotropic  bimodulus  material)  as 
was  shown  in  Chapters  1 1 - IV .  It  is  also  shown  that  the  exact  solutions 
developed  here  offer  a  good  check  for  finite-element  analysis  available 
(see  Chapter  II).  The  results  of  Chapters  III  and  IV  can  also  be  used 
for  this  purpose. 

Good  agreement  was  also  obtained  between  the  approximate 
Galerkin-type  solution  presented  in  Chapter  V  and  an  existing  solution 
for  the  isotropic  bimodulus  case. 

It  has  been  observed  that  for  materials  with  different  proper¬ 
ties  in  tension  and  compression,  the  location  of  the  neutral  surface 
may  vary  considerably  from  the  geometric  mid-plane.  Materials  with 
markedly  different  properties  In  tension  and  compression  are,  of  course, 
most  affected. 

There  is  a  sudden  jump  in  neutral -surface  locations  with  the 
change  In  aspect  ratio  In  the  case  of  free  vibration.  This  is  probably 
due  to  the  eigenvalue  nature  of  this  problem.  They  also  assume  out-of¬ 
plate  values  In  the  case  of  thermal  bending. 

Stacking  sequence  for  two-layer,  cross-ply  laminates  plays  an 
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important  role  In  both  the  amount  of  deflection  and  the  failure  mode. 

For  some  plate  aspect  ratios,  the  maximum  stress  can  be  developed  at  the 
bottom  of  the  plate.  However,  if  the  stacking  sequence  is  reversed,  the 
maximum  stress  may  occur  in  those  fibers  closest  to  the  mid-plane. 

Typical  computer  programs  are  presented  in  Appendix  III. 
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APPENDIX  A 


DERIVATION  OF  THE  PLATE  STIFFNESSES  FOR  TWO-LAYER 
CROSS-PLY  LAMINATE  OF  BIMODULUS  MATERIAL 

In  the  solution  of  problems  involving  laminates  comprised  of 
bimodulus-material  layers,  it  is  necessary  to  evaluate  the  integral  forms 
involved  in  the  definitions  of  the  plate  stiffnesses,  Eq.  (2.1.6).  This 
is  accomplished  here  for  the  case  of  a  two-layer  cross-ply  laminate. 

Each  layer  is  assumed  to  be  of  the  same  thickness,  h/2,  and 
the  same  orthotropic  elastic  properties  with  respect  to  the  fiber  direc¬ 
tion.  Since  each  layer  is  oriented  at  either  0°  or  90°  to  the  x  axis, 
the  laminate  is  also  orthotropic,  i.e.,  there  are  no  stiffnesses  with 
subscripts  16  and  26. 

The  bottom  layer  is  denoted  as  layer  1,  i.e.,  i  =  1  in  Q. 
and  occupies  the  thickness  space  from  z  =  0  to  z  =  h/2,  where  z  is  mea¬ 
sured  positive  downward  from  the  midplane.  The  top  layer  is  denoted  as 
layer  2,  i.e.,  i  -  2,  and  occupies  the  thickness  space  from  z  =  -  h/2 
to  z  *  0. 

In  the  general  case  derived  in  this  derivation,  it  is  assumed 
that  the  upper  portion  of  the  top  layer  (n= 2)  is  in  compression  (k=2  in 
Qijk^)  the  fiber  direction  and  that  the  lower  portion  of  the  top  layer 
is  in  tension  (k=l),  while  the  inner  portion  of  the  bottom  layer  U=l), 
from  z  3  0  to  z  *  znx,  is  in  compression  (k=2),  while  the  outer  portion 
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(from  znx  to  h/2)  of  layer  1  Is  in  tension  (k*l). 

Thus,  the  general  integral  expression  for  A^. ,  the  first  of 
Eqs.  (2.1.6),  may  be  taken  as  the  sum  of  the  integrals  for  each  of  these 
regions: 

Case  1 

Iz^oT  Zx>0) 


rh/2 

V* dz 

rZny 

'  U  QU22  d2  + 


t  nx 


j2l  dz  +  J  Qijll  dz 
znx 


(A-l ) 


Since  the  planar  reduced  stiffnesses  Q-jj^  are  each  respec¬ 
tively  constant  in  the  appropriate  regions,  Eq.  (A-l)  integrates  to  the 
following  result: 

Aij  =  (Qij22  +  Qijll)(h/2)  +  (Qij21  '  Qijll)znx 


+  (Qi j22  "  Qijl2)zny 


(A-2) 


Aij  =  ( 1/2) (Qi j22  +  °ijll^  +  ^1 j21  '  Qijll)Zx 


+  (Q1j22  "  Qijl2)Zy 


Similarly 


(A-3; 


rh/2 

•«  '  J.h/2  112 

f2ny  f0  f2nx  ,h/2 

’  L..  z(,1j22  dz  *  J  ;(!1J12  dz  *  J  2<51j2)  dz  *  J  zQtjll  dZ 


(A-4) 


'  *  <W(h2/8)  *  «1j2,  '  Q1jll)<2nx2/2) 

*  ((!U22  ’  Q1jl2>  V/2) 

or 

BU/h2  •  0/8)(-  «1J22  *  Qun>  *  <1fJ21  •  QtjM»(V/2> 

*  «IJ22  -  «ij,2><y'2) 

Also 


°ij 


rh/2 

J-h/2 


2%l«  dz 


ny 

f° 

fnX 

z2Q1 122  dz  + 

-h/2 

L  2%12d2t 

l„  z2Q1j21  dl*. 

ny 


h/2 


z2Qijll  dz 


nx 


*  (Q(j22  +  Qi j, l><h3/24>  MQ(J21  -  Q,j„)(V/3) 

*  (Q1 J22  '  Qijl2)(z„y3/3> 

or 


(A-5) 


(A-6) 


(A-7) 


(A-8) 


Dfj/h3  .  (l/24)(Q1j22  *  qtj„,  *  (QiJ2i  -Q,j„)(ZxV3) 

(A-9) 

+  (Qij22  ’  Q1j12^Zy3/3) 

Similarly 

Case  2 

(zy>0.  zx<0) 

rh/2 

Au  ‘  l.w  Quu  82 

f2™  rO  ,zny  .h/2 

'  U/2  5,322  iZ  *  It  5,3,2  di  *  Jo  Q,J2'  iZ  *  J,  Q1JU  82 

nx  *ny 
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=  (Qijll  +  Qij22)(h/2)  +  (Q1j22  ‘  Qijl2)znx  +  (Qij21  "  Qijll)zny 


or 


A1j/h  =  (Qijll  +  Qij22)/2  +  (Qij22  ‘  Qijl2)Zx  +  (Qij21  '  Qijn)Z} 

and 


(A— 10) 


BU/h2  ■  (Q1jtI  -  QiJ22)/8  +  «1J22  -  Q1jl2>‘Zx2/2>  +  «1J21  - 

V’  '  (Q1jll  *  0,j22)/24  ♦  (QjJ22  -  Qjjl2)(2xV3)  *  (QfJ2|  -  Q,j„)(Z//3) 


Case  3 

TVU7  zy>0) 

V  ’  +  Q«22)/2  +  (QU21  -  Qijll)Zx 

Bi/h2  ■  «1J11  •  Bij22^ZB  *  «1j21  -  <!fjll)<Zx2'2>  (A-m 

VM  ’  (,1J11  *  W/24  +  (Q1J21  '  W(V/3) 

Case  4 

TT^T  zy<o) 

A1j/h  ‘  «ijU  *  B1j22>Z2  +  <Bij22  -  8ijl2>Zy 

B,j/ha  ■  <Q(ji]  "  Bij22>/8  *  <Blj22  "  Bljl2*(Zy2/2>  (A-12) 

V”3  .  (Qu1]  *  QjJ22)/24  ♦  (QfJ22  -  Q,j,2)(2y3/3) 


for  neutral  surface  going  out  of  plane. 


For  a  single  layer  (0°), 


A1j/h  "  (Qijll  +  Qij21)/2  +  (Qij21  “  Qijll)Zx 
B1j/h2  *  (Q1j11  ’  Qij21)/8  +  (Q1j21  "  Q1jll)Zx/2 
°1j/h3  *  (Q1jll  +  Qij21)/24  +  (Q1j21  ’  Q1jll)Zx/3 
For  neutral  surface  out  of  plane. 


Case  1 

TT^TS) 


V  '  Qij21 

V"2  ■  0 
VM  ■  QijZl/12 


Case  2 

TT^.5) 


A1j/h  "  Q1jll 
B1j/h2  "  0 

Vh3  *  Qijn/12 


(A- 17) 


(A-18) 


(A- 19) 
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APPENDIX  B 

DENSITIES  OF  ARAMID,  POLYESTER,  AND  RUBBER 


Fibers: 

According  to  the  Kevlar  Data  Manual  [50],  the  densities  of  the 
fibers  of  aramld  and  polyester  are: 

Kevlar  (aramid)  *  0.052  Ib^/in3 

Polyester  ■  0.049  lb^/ln3 

Rubber: 

Reference  [51]  lists  nautral  rubber  and  isoprene  rubber  at 
specific  gravity,  0.93,  and  SBR  plus  BR  rubber  at  specific  gravity,  0.94. 
Taking  0.93,  the  density  of  rubber  can  be  calculated  as: 

Density  «  0.93  x  62.4  lbf/ft3  (H2Q)/1728  in3/ft3 

■  0.034  lbf/in3 

Composite: 

If  p  represents  the  density  (f  for  fiber  and  m  for  matrix)  and 


V  represents  the  volume  fraction,  the  density  of  the  composite  p  an  be 
written  as: 


The  values  of  Vf  according  to  Ref.  [2]  are: 

Aramid/ rubber  :  0.140 
Pol yes ter/ rubber:  0.149 

Thus,  the  densities  of  the  composites  can  be  listed  as  follows 


Composite 

Arami d/rubber 
Polyester/ rubber 


Density 

p(lbf/1n3) 

0.037 

0.036 


APPENDIX  C 

THERMAL  FORCE  AND  MOMENT  EXPRESSIONS  FOR  CASES  II  -VIII 
T  T  T  T 

Expressions  for  N1|X,  N2<y,  M1>x,  and  M2>y  were  derived  for  Case 
I  in  Chapter  IV.  In  a  similar  way,  one  can  obtain  the  expressions  for  the 
above-mentioned  quantities  for  the  remaining  seven  cases  as  follows: 


Case  II 


z  >0, 
nx  ’ 


ny 


>0 


Nllx  *  at^122  +  6111)^0^/2)  +  (®121  *  8m)^oznx^ 

+  (Sin  -  6i22)(T1h/8)  +  (Bi2i  -Bin)(TiZn2x/2h)> 

N2^y  *  B{(B222  +  B2n)(TQh/2)  +  {8221  '  62n)(Toznx) 

+  (8211-8222)^1^/3)  +  (8221  ■  82ii)(^i*n2x/2h)} 

(C-l) 

Milx  *  “UBih  -  8122)(T0h2/8)  +  (8121  "  6m)(t0Zn2x/2) 

+  (8122  +  8m)(T1h2/24)  +  (8121  ■  BmHTjZ^/Sh)} 

M2,y  »  81(8211  “  8222 ) (^q^z/3)  +  (6221  "  621 1 ) (^0zn2x/2) 

+  {8222  “  62 1 1 ) ( ^ 1 h2/24)  +  {$221  “  62ii)(Tj z  3 / 3h ) > 
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Ni.x  ”  at(Si22+6n i ) (TQh/2)  +  (8121  -  8111 )(^0zny) 

+  (8122  ■  Sii2)(f0Znx)  +  (^m  ■  6122 )  (^  lh/8) 

+  (8121  -  8111  )(fizny/2h)  +  (8122  “  8n2)(TiZn2x/2h)} 
N2fy  *  8{(8222  +  ^21l)(^Q*,'/2)  +  (8221  ■  S21l)(^0zny) 

+  (8222  "  ®212)(^0Znx)  +  (8211  ~  8222)^1^/8) 

+  (0221  ■  8211  )(flZn2y/2h)  +  (8222  ■  8212)(^lzn2x/2h)  ) 
Mllx  *  a{(8iu  -  6i22)^0hz/8)  +  (8121  '  8m )(TQ2n^/2) 

+  (8122  "  8n2)(T0zn2x/2)  +  (8122  +  8m)(T1h2/24) 

+  (8121  ■  8ni)(fizn3y/3h)  +  (8122  ’  8n2)(^izn3x/3h)} 

M2,y  *  8{(82n  -  0222)^oh2/8)  +  ( B221  "  82 1 1 ) (^Z^/2) 

+  (8222  -  62i2)(Vnx/2)  +  (8222  +  82i  1 )  (Tj  h2/24 ) 

+  (8221  "  82li)(Ti2n^/3h)  +  (8222  “  82l2)(TiZn3x/3h)} 


Case  IV 

z  <d.  z  <0 
nx  *  ny 

Nijx  *  ®t (8122  +  Sin) (TQh/2)  +  (8122  -  Si2i)(T0zny) 

♦  (8111  -  0 1-22 ) (^ l^*/8)  +  (8122  *  8i2i)(TiZn:y/2h)} 
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N^y  *  6t(&222  +  ^2 1 1 ) (^0^/2)  +  (6222  -  &22l)(^0Zny) 

+  (8211  -  ^222 ) ( ^1 ^/8 )  +  (0222  "  ^221 ) (f 1 zny/2/h) } 

Mjx  3  ( 6 1 1 1  ■  S122 ) (TQh2/8)  +  (8122  “  ^122 )(^0zny/2) 

+  (8^22 +  8m)(T^h2/24)  +  (6  122  "  &mKfiZn3y/3h)} 


M2^y  *  8((0211  ■  8222)(Tgh2/8)  +  (8222  •  6221 )(^0zny/2) 

+  ($222  +  S21l)(Tj h2/24 )  +  (8222  “  ^221 ) ^lzny/^^)^ 
For  neutral  surface  going  out  of  plane, 

Case  V 


I7O,  Zy<-0.5 

Njx  a  aH0ui  +  ^ii2)^0/2  +  (8121-8112)^1/8} 

N2,y  *  8( (8221  +  8212)^q/2  +  (8221-8212)^1/8} 
Mi^x  3  a{ (0121  "  8112)^/8  +  (8121  +  13 1 1 2 ) ^ i/24 } 
^zly  3  01(8221  ’  8212)^5/®  +  ^®221  +  8212)^1/24} 


Nilx  3  «<(0121  ♦  8112)^/2  +  (8121-8112)^1/8} 

N2^y  3  8(  (0221  +  8212)^/2  +  (8221  -  8212)^1/8} 

MiTx  “  <>{(8121  -  8112)^/®  +  (8121  +  8h2)Ti/24} 

Mj^y  3  8{(8221  -  8212)^/®  *  (8221  +  8212)^1/24} 


(C-3) 


(C-4) 


(C-5) 
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Case  VII 
?x>0.5,  ly>Q.S 

"ilx  *  a{^ni  +  @h2)t0/2  +  (sln  -  g112)Ti/8} 
Nz,y  *  8{(e2n  +  82i2)Tq/2  +  (e2u  -  8212)Tl/8> 
"ilx  *  a{(8m  -  8in)V8  +  (0ni  +  Si12)T1/24} 
M2,y  *  e{(e211  -  8212)To/8  +  (e211  +  e212)fj/24} 


Case  VIII 

^0.5,'Z. 


<-0.5 


Nl«x  *  a^8ul +  8U2)To/2  +  (8ln  -  8U2)Tx/8} 
N2ly  =*  8{(8221  +  8222)To/2  +  (6221  -  8222)T1/8} 
Ml!x  *  a{(s121-6l22)T0/8  +  (8121  +  8122)V24} 
M2,y  *  6{(8221  -  8222)To/8  +  (8221  +  8222)Tj/24} 


(C-6) 


(C-7) 


For  a  single  layer,  change  a112  to  8ln,  e122  to  812I,  s212  to  s2u 
and  8222  to  822i. 
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APPENDIX  D 

IN-PLANE  DISPLACEMENTS  AND  SLOPE  COEFFICIENTS 

The  in-plane  displacements  and  the  slope  coefficients  correspond¬ 
ing  to  the  middle-surface  deflections  presented  in  the  preceding  chapters 
are  listed  in  the  following  tables. 


Table  D2.5.  In-Plane  Displacements  and  Slope  Coefficients 
Associated  with  Table  2.5 


a/b  (UE22ch3/qob4)xl0"2 

(VE22' 

ch3/qob4)xl0'2 

(XE22( 

:h3/qob4)xl0'2 

(YE22‘ 

:h3/qob4)xl0'2 

Aramid-Rubber  (NL=1) 

0.5 

0.0510 

• 

0.0238 

- 

0.1162 

- 

0.0730 

0.6 

0.0824 

- 

0.0386 

- 

0.1883 

- 

0.1333 

0.7 

0.1204 

- 

0.0552 

- 

0.2754 

- 

0.2187 

0.8 

0.1725 

- 

0.0720 

- 

0.3727 

• 

0.3295 

0.9 

0.2064 

- 

0.0877 

- 

0.4745 

- 

0.4634 

1.0 

0.2495 

- 

0.1013 

- 

0.5754 

- 

0.6161 

1.2 

0.3264 

- 

0.1209 

- 

0.7575 

- 

0.9569 

1.4 

0.3844 

- 

0.1308 

- 

0.8982 

- 

1.3101 

1.6 

0.4222 

- 

0,1313 

- 

1.0509 

- 

1.9486 

1.8 

0.4427 

- 

0.1269 

- 

1.0778 

- 

2,2132 

2.0 

0.4503 

- 

0.1269 

- 

1 .0778 

- 

2.2132 

Polyester-Rubber  (NL»1) 

0.5 

0.0150 

- 

0.0073 

- 

0.0503 

- 

0.0449 

0.6 

0.0252 

- 

0.0122 

- 

0.0848 

- 

0.0787 

0.7 

0.0387 

- 

0.0184 

0.1304 

- 

0.1281 

0.8 

0.0554 

- 

0.0256 

- 

0.1871 

- 

0.1960 

0.9 

0.0749 

- 

0.0335 

- 

0.2539 

- 

0.2845 

1.0 

0.0967 

- 

0.0419 

- 

0.3291 

- 

0.3945 

1.2 

0.1442 

- 

0.0588 

- 

0.4947 

- 

0.6760 

1.4 

0.1915 

- 

0.0744 

- 

0.6616 

- 

1.0211 

1.6 

0.2327 

- 

0.0873 

- 

0.8096 

- 

1.3972 

1.8 

0.2646 

- 

0.0971 

• 

0.9264 

- 

1.7716 

2.0 

0.2863 

- 

0.1041 

- 

1.0087 

- 

2.1199 

Table  D2.6.  In-Plane  Displacements  and  Slope  Coefficients 
Associated  with  Table  2.6 


a/b  (UE22ch3/qob4)xl0'2 

(VE22 

ch3/qob4)xl0'3 

(xe22( 

ch3/qQb4)xl0'2 

(ye22 

ch3/qob4)xl0‘2 

Aramid-Rubber  (NL=2) 

0.5 

0.5967 

- 

0.5570 

- 

0.1114 

- 

0.0716 

0.6 

0.4608 

- 

0.4168 

- 

0.1789 

- 

0.1290 

0.7 

0.3612 

- 

0.3197 

- 

0.2602 

- 

0.2099 

0.8 

0.2855 

- 

0.2487 

- 

0.3517 

- 

0.3151 

0.9 

0.2267 

- 

0.1951 

- 

0.4485 

- 

0.4429 

1.0 

0.1804 

- 

0.1539 

- 

0.5456 

- 

0.5899 

1.0 

0.1804 

- 

0.1539 

- 

0.5456 

- 

0.5899 

1.2 

0.1148 

- 

0.0969 

- 

0.7247 

- 

0.9220 

1.4 

0.0736 

- 

0.0618 

- 

0.8672 

- 

1.2718 

1.6 

0.0478 

- 

0.0400 

- 

0.9673 

- 

1.6084 

1.8 

0.0315 

- 

0.0265 

- 

1.0289 

- 

1.9144 

2.0 

0.0211 

- 

0.0179 

- 

1.0603 

- 

2.1832 

Polyester-Rubber 

(NL=2) 

0.5 

0.5139 

- 

1.1570 

- 

0.1169 

- 

0.0739 

0.6 

0.3980 

- 

0.9203 

- 

0.1881 

- 

0.1330 

0.7 

0.3113 

- 

0.7477 

- 

0.2731 

- 

0.2152 

0.8 

0.2440 

- 

0.6133 

- 

0.3660 

- 

0.3195 

0.9 

0.1909 

- 

0.5045 

- 

0.4597 

- 

0.4417 

1.0 

0.1489 

- 

0.4146 

- 

0.5475 

- 

0.5754 

1.2 

0.0898 

- 

0.2787 

- 

0.6882 

- 

0.8497 

1.4 

0.0541 

- 

0.1871 

- 

0.7730 

- 

1.0991 

1.6 

0.0330 

- 

0.1267 

- 

0.8099 

- 

1.3047 

1.8 

0.0206 

- 

0.0872 

- 

0.8138 

- 

1.4662 

2.0 

0.0132 

- 

0.0612 

- 

0.7979 

- 

1.5905 

Table  04.2.  In-Plane  Displacements  and  Slope  Coefficients 
Associated  with  Table  4.2 


a/b  Uh/c^  Ttb2  (Vh/^1  Tib2)xl01  Xh/aa*  Txb2  (Yh/cta*  T1b2)xl01 


Aramid-Rubber  (NL=1) 

0.5 

- 

0.1596 

-  0.1163 

-  6.9363 

-  0.6982 

0.7 

- 

0.2226 

-  0.1669 

-  0.1430 

-  0.6081 

0.9 

- 

0.2854 

-  0.2034 

-  0.2395 

-  0.5857 

1.0 

- 

0.3168 

-  0.2172 

-  0.2959 

-  0.5886 

1.2 

- 

0.3795 

-  0.2383 

-  0.4236 

-  0.6131 

1.4 

- 

0.4421 

-  0.2531 

-  0.5692 

-  0.6556 

1.6 

- 

0.5045 

-  0.2638 

-  0.7298 

-  0.7110 

1.8 

- 

0.5669 

-  0.2716 

-  0.9018 

-  0.7765 

2.0 

- 

0.6292 

-  0.2775 

-  1.0813 

-  0.8495 

Polyester- 

•Rubber  (NL=1) 

0.5 

- 

0.1619 

-  0.1333 

-  0.0767 

-  0.2048 

0.7 

- 

0.2226 

-  0.1848 

-  0.1587 

-  0.1833 

0.9 

- 

0.2826 

-  0.2197 

-  0.2639 

-  0.1846 

1.0 

- 

0.3124 

-  0.2325 

-  0.3238 

-  0.1906 

1.2 

- 

0.3719 

-  0.2516 

-  0.4550 

-  0.2108 

1.4 

- 

0.4310 

-  0.2647 

-  0.5959 

-  0.2393 

1.6 

- 

0.4898 

-  0.2740 

-  0.7388 

-  0.2736 

1.8 

- 

0.5482 

-  0.2807 

-  0.8760 

-  0.3114 

2.0 

- 

0.6062 

-  0.2858 

-  1.0002 

-  0.3503 
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Table  04.3.  In-Plane  Displacements  and  Slope  Coefficients 
Associated  with  Table  4.3 


a/b 

(Uh/a2t  T1b2)xl0'1 

(Vh/a2t  T1b2)xl0"1 

Xh/a2t  T^2 

Yh/a2t  Ttb2 

Aramid-Rubber  (NL*2) 

0.5 

-  3.3455 

-  0.1701  x  104 

2.6210 

-  1.1264 

0.7 

-  5.4763 

-  0.2371  x  104 

4.7012 

-  1.0100 

0.9 

0.0229 

-  0.7557 

-  0.0397 

-  1.2135 

1.0 

0.0547 

-  0.7789 

-  0.0558 

-  0.2227 

1.2 

0.1112 

-  0.8233 

-  0.0853 

-  0.2403 

1.4 

0.1513 

-  0.8597 

-  0.1081 

-  0.2548 

1.6 

0.1738 

-  0.8867 

-  0.1240 

-  0.2655 

1.8 

0.1821 

-  0.9054 

-  0.1341 

-  0.2730 

2.0 

0.1803 

-  0.9177 

-  0.1401 

-  0.2779 

Polyester-Rubber  (NL=2) 


0.5 

-  46.395 

- 

32.850 

6.2344 

-  3.9037 

0.7 

-  66.708 

- 

33.152 

8.9231 

-  3.8930 

0.9 

-  63.472 

- 

34.827 

8.1669 

-  4.2283 

1.0 

-  55.893 

- 

36.025 

6.9217 

-  4.5183 

1.2 

-  40.106 

- 

37.695 

4.3247 

-  4.942 

1.4 

-  37.079 

- 

0.1555 

X 

103 

1.1373 

0.3559 

X 

103 

1.6 

-  42.873 

- 

0.1650 

X 

103 

1.2790 

0.3544 

X 

103 

1.8 

-  47.826 

- 

0.1578 

X 

103 

1.4415 

0.3404 

X 

103 

2.0 

-  52.199 

- 

0.1428 

X 

103 

1.6249 

0.3218 

X 

103 

Table  05.2.  In-Plane  Displacements  Associated  with  Table  5.2 


Aramid-Rubber  (Nl«l 


Polyester-Rubber  (NL*1) 


0.5  0.7656  x 
0.6  0.1195  x 
0.7  0.1605  x 
0.8  0.1857  x 
0.9  0.1816  x 
1.0  0.1382  x 
1.2  -0.6104  x 
1.4  -0.3813  x 
1.6  -0.7777  x 
1.8  -0.1216  x 
2.0  -0.1675  x 


-0.3512  x 
-0.6279  x 
-0.1012  x 
-0.1543  x 
-0.2295  x 
-0.3351  x 
-0.6757  x 
-0.1221  x 
-0.1981  x 
-0.2949  x 
-0.4106  x 


0.2051  x  10 
0.2072  x  10‘2 
0.9315  x  10'3 
-0.1215  x  10 
-0.4870  x  10 
-0.9015  x  10 
-0.1856  x  10 
-0.2876  x  10 
-0.3886  x  10 
-0,4835  x  10 
-0.5686  x  10‘ 


-0.8036  x 
-0.1289  x 
-0.2085  x 
-0.3786  x 
-0.5964  x 
-0.9549  x 
-0.2047  x 
-0.3620  x 
-0.5615  x 
-0.7949  x 
-0.1053 


Table  05.3.  In-Plane  Displacements  Associated  with  Table  5.3 


0.5 

0.4194 

X 

10"5 

-0.6301 

X 

10-6 

0.1982 

X 

10'4 

-0.9016 

X 

10"5 

0.6 

0.7255 

X 

10'5 

-0.1140 

X 

10'5 

0.3383 

X 

10"4 

-0.1823 

X 

10"4 

0.7 

0.1148 

X 

10'4 

-0.1889 

X 

10'5 

0.5280 

X 

10'4 

-0.3298 

X 

10"4 

0.8 

0.1705 

X 

10‘4 

-0.2944 

X 

10"5 

0.7696 

X 

10“4 

-0.5484 

X 

10"4 

0.9 

0.2412 

X 

10'4 

-0.4382 

X 

10"5 

0.1062 

X 

10'3 

-0.8536 

X 

10‘4 

1.0 

0.3282 

X 

10"4 

-0.6299 

X 

10"5 

0.1399 

X 

10'3 

-0.1259 

X 

10'3 

1.2 

0.5547 

X 

10"4 

-0.1200 

X 

10"4 

0.2173 

X 

10"3 

-0.2402 

X 

10‘3 

1.4 

0.8529 

X 

10'4 

-0.2126 

X 

10‘4 

0.2972 

X 

10'3 

-0.3993 

X 

10‘3 

1.6 

0.1214 

X 

10’3 

-0.3567 

X 

10"4 

0.3708 

X 

10"3 

-0.5964 

X 

10'3 

1.8 

0.1626 

X 

10“3 

-0.5737 

X 

10"4 

0.4308 

X 

10‘3 

-0.8177 

X 

10‘3 

2.0 

0.2055 

X 

10'3 

-0.9012 

X 

10"4 

0.4750 

X 

10'3 

-0.1046 

X 

10~2 

0.5  0.2483  x  10’2  -0.1847  x  10“2  0.7623  x  10"3  -0.1374  x  10‘3 
0.6  0.2562  x  10"2  -0.3998  x  10“2  0.1194  x  10'2  -0.2888  x  10"3 
0.7  0.1401  x  10"2  -0.7835  x  10'2  0.1599  x  10"2  -0.6180  x  10"3 
0.8  -0.1053  x  10"2  -0.1391  x  10"1  0.1849  x  10"2  -0.1338  x  10'2 
0.9  -0.4506  x  10'2  -0.2240  x  10”1  0.1807  x  10'2  -0.2729  x  10'2 
1.0  -0.8479  x  10“2  -0.3307  x  10_1  0.1391  x  10‘2  -0.5014  x  10"2 
1.2  -0.1568  x  10’1  -0.5676  x  10”1  -0.4823  x  10"3  -0.1241  x  10_1 
1.4  -0.1938  x  10_1  -0.7625  x  10'1  -0.2872  x  10“2  -0.2195  x  10_1 
1.6  -0.1973  x  10” J  -0.8835  x  10"1  -0.2872  x  10’2  -0.2195  x  10'1 
1.8  -0.1835  x  10'1  -0.9482  x  10'1  -0.5559  x  10"2  -0.3967  x  10_1 
2.0  -0.1636  x  10"1  -0.9814  x  10"1  -0.5559  x  10'2  -0.3967  x  10'1 
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APPENDIX  E 

TYPICAL  COMPUTER  SOLUTIONS 

Sample  computer  printouts  are  attached.  The  programs  are 
presented  separately  for  each  problem.  They  are  self-explanatory.  The 
dimension  statement  allows  required  values  for  different  quantities.  For 
example,  reduced  stiffnesses  take  5x5x2x2  values  (2  represents  2  layers, 
and  the  other  2  Indicate  different  properties  In  tension  (1)  and  compres¬ 
sion  (2))  and  THETA  takes  either  zero  or  1.5708  radians. 

The  ISML  library  subroutine  LEQT2F  solves  N  equations  In  N 
unknowns.  EIGZF  solves  eigenvalue  problem,  and  ZRPOLY  solves  cubic 
equation.  It  should  be  mentioned  that  WKAREA,  WK,  etc.  are  part  of  such 
subroutines. 

The  values  of  z„w  and  z_  are  assumed  In  the  beginning  to  get 
nx  ny 

the  displacements.  An  Iterative  procedure  Is  then  adopted  until  a  pre¬ 
cision  of  +  0.0001  Is  achieved  to  get  the  actual  deflections.  Later  on 
the  whole  procedure  Is  repeated  for  different  aspect  ratios. 

Computations  are  carried  out  for  a  single-layer  rectangular 
plate  by  making  NL-1  (number  of  layers)  and,  of  course,  with  a  proper  set 
of  equations  for  A^,  B^,  and  0^.  With  proper  A^,  B^,  and  and 
NL«2,  the  program  calculates  the  deflections  and  neutral -surface  loca¬ 
tions  for  a  two-layer,  cross-ply  rectangular  plate.  THETA  takes  the 
values  0  and  1.5708  (corresponding  to  0°  and  90°)  radians  for  the  two- 
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COMPUTER  PROGRAM  FOR  STATIC  BENDING  (SINGLE  LAYER) 


»  HBXWWWW H«» -  xmr+Kj***.* 


TItfE=( 1.5E) 

IMPLICIT  Sc  JL«  :j<  A-r-  .  C  ~L  ) 

0  ?  ■  E  f-  S  ic.\  Z  1  J  (  2  )  .  0  2  t  l  2  )  *  jl  •;{  2  #  iC?.‘i(C  )  .OJi  !?i.  4Nvl2  <  2  )  .  »  MJ2  1(2)  *  : 

*G(5.S»2}.C{£,E.:>.2).-«Z.  ;■  )  ,«.•  (  j  ,  1  >  ,0  (  ,  ,  j  ,<  t  3  ,  f, )  ,  CD  A  EH  t  t  .  1  i  , 

P t=  >i  (€0)  .  Ti-cTA  !2  i 

SE.-C.2  1  H  n.:“ll(2!.i2c{U»  222 i 2)  .  i>  1  2  {  l)*Gl?{2)  »  w  ’  j  {  I  }  •  6  C  2  ( 2  )  • 

'  u  3  1  (  *  /  if  1  j  i  1  ^  2  J  •  2itl)*ANv3t*6f#  T  i~  c.  T  A  i  l  l  1 

*lr:£TAU) 

C  C  2D  1=1.5 
o  c  ;:•  o  j  a  i.s 
C  2  G  K  -  1 .  t 
CC  20  L -  1  *  2 
GQ ( I  .  J ,K  )  =  0 . 3 
0<i«J.K,L>*0.o' 

1  CCfTIMjH 
N  L  —  1.0 

30  1=  1  .  £ 

CC  20  „=  1  » £ 

CC  30  K=  1, 1 

OEf^CV  a  1-ANLlt  (  K  JtAMIt  1{K  » 

(ia(t.i.K)sen(K)/csNCN 

DC  (2  •  2  .  K  )-e<  2  (K  )  /CENC?' 

CQ(l,2,K)s4M,12(X)tCi.  (2.2  ,  K ) 

GC  (  1  .K  )  -OC  (  1, 2  »K  ) 

CO  <3 .3  ,K  )  =  Cl  2  iK  ) 

<10(4 ,  4  ,K  >  =  C2  3(K) 

CC  <S  *  2  .  K  )*C2  l  <*.  > 

DC  30  L  =  l  ,  M. 
x=CCCS(  Tr-CTA  (L  )  ) 

Y  =  CS!MTHSTML)> 

X2  =  K*  ■>  2 

Y  2=>  *  42 
X  4  =  X  «•  4  4 
Y4  = 

G  (  l  *  1>  K  .  L  J  =  CD  '  t  .  1  .  K  )  »  X  4*  CO  (  2 . 2  »  K  J  4  A  4 

C(.l.2.K,L)  =  C  0  (  I  «  2  /  K  )  <•  (  Y  4  >  X  4  > 

C(2»1.K,L)=C<  1  *  2 . K  ,  .  ) 

C  (2  »2  >K  it  )  =  C Q  (  l  «  1  ,  K  1  4  Y <t  4  CO  (  <2  *  2  « <0  *  >  *t 

JlwtK  I(<)  <  (  V  4  1  A  A  < 

0(4»<t»Kil)sC5(4«4*K'(x?iCiD(S*5«K14’2 

C  (S  .5  .K  .  L  r-OC  (4 ,4  ,  S  )<i  Y2  +  CG  ij  ,  j  ,  ><  )*>.2 

CCNT INLc 

Ax=5.0 

a  r=  t  o  .o 

2  X*0 •  OC  1 

ZY*-C .00 i 

PRINT  16.2tl.E22.G12.G22.Gol  .  4  NU  l  2  .  A  >  C2  1  .  TliET  4  .Act*  .T3AH0  »  TEAS!. 
Zx2  =  Zx»  * 2 
2,:2-Z't»4  2 
2  >  2s  2  ;<  •*  *  2 
Z Y3*ZY  »*2 
wvs  l 
NN*1 

alpha  a  pv43#1416/ax 
OCT  4  a  NN<-  0  .  1*1  16/£T 
CC  5  C  1=1,2 
OC  3  C  2=1.2 

A  <  1  ,  J  1  =  |C(  I...  i.  I)  m;(  I  JiZi  1  )  )/?HC(li  J.Si  1  !*,:H  !  .Jil  .  1  )  >  *  l  X 

o<  i .  j  i  *  to  ( i ,  j,  i .  i )  -<;  ( i .  j  .2,  i  >  i/attu  i.  j.  c .  i  >  -<;<  i.j.i.n>*zx2/s 

0(  1  .  J  >*  (  0  (  1  .  >  .  t  .  1  )  40  (  l  .2  .  n  1  /24M  0  (  l  .  J.2 . 1  >  -  C(  I  .  J  .  1  .  I  )  >*ZX  J/3 


»  voV-Si.' 


21 -•  IUMI 
C(j.IMCU.J) 

CCM  INCH 

» <• «  *  i o  i 4  . <» .  i .  n  ♦  <  •» ,  i  , c .  ’ '  :  i  i  )••  1. 1  * ;  .  i  )  >  .; « ♦ 

>  ( c  ( *  • »» »  z. .  3  >  - : i a  •  o  •  i ,  ?  ; >  -  *  f 
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‘A  differential -equation  formulation  Is  presented  for  the  equations  governing 
the  small -deflection  elastic  behavior  of  thick  plates  laminated  of  anisotropic 
bimodulus  materials  (which  have  different  elastic  stiffnesses  depending  upon 
the  sign  of  the  fiber-direction  strains).  This  study  Is  concerned  with  the 
problems  of  bendlno  due  to  pressure  loading,  free  vibrations,  and  thermal  ex-  ■ 
pension  of  thick  plates  having  finite  transverse  shear  moduli.  Large  static 
deflections  of  thin  plates  are  also  analyzed  bv  an  approximate  technique.  A  , 
simple  scheme  Is  used  for  determining  explicitly  the  location  of  the  (over) 
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20.  Abstract  -  Cont'd 

^neutral  surface  which  Is  defined  on  the  basis  of  fiber-direction  strain.  This 
provides  a  basis  for  rational  design  of  such  plates.  Exact  closed-form 
solutions  are  presented  for  freely  supported  rectangular  plates  subjected  to  a 
sinusoidally  distributed  normal  pressure.  Based  on  experimentally  measured 
bimodulus  properties,  some  numerical  computations  are  carried  out.  Good 
agreement  Is  obtained  when  compared  with  numerical  results  existing  In  the 
literature  for  special  cases.  The  results  presented  here  can  be  used  to 
validate  finite-element  codes  being  developed  for  analysis  of  thick  plates 
laminated  of  bimodulus  material.  «=• 
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